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ABSTRACT 


Various properties of multivariable root loci are analyzec! from a 
frequency-domain point of view by using the technique of Newton polygons, 
and some generalizations of the SISO root loctis rules to the multivariable 
case are pointed out. The behavior of the angles of arrival and departure 
is related to the Smith-MacMillan form of G(s), and explicit equations for 
these angles are obtained. After specializing to first-order and a restricted 
class of higher-order poles and zeros, some simple equations for these angles 
that are direct generalizations of the SISO equations are found. 

The unusual behavior of root loci on the real axis at branch points is 
studied. The SISO root locus rules for break-in and break-out points are 
shown to generalize directly to the multivariable case. Some methods for 
computing both types of points are presented. 

An equation for the number of loci on the real axis at any point is 
derived. The special cases of 2)<2 G(s) and symmetric G(s) are investigated 
separately. Finally, for high gains, equations for the first-order asymptotes 
and pivots are derived, and previous results on higher-order asymptotes are 
reviewed. 
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CHAPTER I 
Introduction 


1.1 Ntotivation and Sunaiary of Results 

The main aspect of the problem of designing a feedback compensator 
for a linear system is that one seeks to characterize the effect of 
feedback on the closed-loop behavior of the system. A balance must be 
struck between a characterization so simple that no real insight into 
the closed-loop system behavior is gained, and one so complex that 
interpretation of it is too difficult to be viseful. 

The root locus does a reasonable job of striking such a balance. The 
root locus teclinique consists of plotting the paths of the movements of 
the closed- loop poles in the complex plane as a single feedback parameter 
is varied. This has two advantages: the locations of the closed- loop 

poles furnish considerable information on the response of the system, 
particularly the transient response; and variation of a single parameter 
gives crude, but simple, notions of the options in pole assignment and 
the identities of specific poles that must be shifted. 

The major disadvantage of the root locus is that little information 
is furnished on robustness of the closed- loop system. Thus the root locus 
method nicely complements stability tests such as the Nyquist stability 
criterion, which give information on system robustness (e.g. phase and 
gain margins) but only "yos-or-no" information on stability and no infor- 
mation at all on the form of system responses. It should also be noted 
that by suitable reformulation the effect of variation of an uncertain 
system parameter can be studied using the root locus. 


The root locus for single- input -single-output (SISO) systems was 
first studied by Evans in 1948 [14] and is treated in detail in any 
decent elementary control theory text (e.g. [7], [8]]: a set of rules 
for plotting it is given in Chapter II. Investigation of the root locus 
for multivariable systems began in earnest in the mid-1970s, and is far 
from being concluded. Research has proceeded along two main lines: the 

state-space approach used by Shaked, Kouvaritakis , and Owens (e.g. [9], 
[IS], [16]) and the frequency-domain approach used by Postlethwaite and 
MacFarlane (e.g. [1], [2], [3]). 

The state- space approach seems to have been better suited for inves- 
tigating the behavior of root loci that approach infinity for high feed- 
back gains, and a considerable body of knowledge has been amassed on this 
subject. However, less information is available on the angles of arrival 
and departure of loci at zeros and poles, and almost none on breakin and 
breakout points and the presence of loci on portions of the real axis. 

The pioneering work of Postlethwaite and MacFarlane suggests that their 
frequency-domain methods may be better suited for investigating these 
issues. 

The aims of this thesis are threefold: first, to develop and extend 

the frequency-domain methods of Postlethwaite and MacFarlane into explicit 
results and equations for the angles of arrival and departure and for the 
locations of loci on the real axis; second, to show how the well-known 
SISO root locus rules do or do not generalize to the multivariable case; 
and third, to serve as a compendium of rules for plotting the multivariable 
root locus. 

We start off by laying in Chapter II the groundwork for the material 


to follow. The basic problem is presented, and basic equations obtained 
from it. Various features of the root locus, such as poles, zeros, 
angles, and Butterworth patterns, are defined, and the SISO root locus 
rules are reviewed. 

We will obtain results on angles of loci by finding series approxi* 
mat ions for the loci near the point of interest. These will be obtained 
by using the Newton polygon technique, an ingenious graphical device 
presented and demonstrated in Chapter II. 

Results for angles of arrival and departure are presented in Chapter 
III. After quickly reviewing the state-space results on this subject, new 
results and explicit equations for the cases of both first-order and 
higher-order poles and zeros are obtained. In particular, the case of 
higher-order poles and zeros turns out to be vastly more complex than 
might be expected from the SISO rules. 

After considering the general case, we specialize to first-order 
and a certain ’’simple" class of higher-order poles and zeros. Using an 
approach based on constructing a Laurent series of the system transfer 
function, simpler equations are obtained that turn out to be nice 
generalizations of the SISO equations. These results have not, to our 
knowledge, appeared in the literature. 

Unusual behavior rf the root locus, such as a locus on the real axis 
suddenly turning around, are associated with entities called branch 
points. We discuss these briefly in Chapter IV, and present some equations 
describing their effects on root loci. Also in Chapter IV, breakin and 
breakout points of multivariable root loci are investigated for the 
first time. Equations are obtained for computing them, and the angles of 


10 


loci breaking in or out are shown to be the same as in the SISO case. 

In Chapter V« we address the previously uninvestigated problem of 
determining the location of root loci on the real axis. This turns out 
to be vastly more complex than the SISO case, since more than one branch 
can lie on the real axis at a given point. An equati^ is obtained for 
the general case, and other results presented. The case of two inputs 
and two outputs is investigated separately, and is shown to be considerably 
simpler than the general case. 

In Chapter VI the asymptotic behavior of root loci is considered. 

For the generic case of first-order asymptotes and pivots, simple 
equations are given for angles and pivots. For higher-order asymptotes, 
the results of Kouvaritakis and Shaked [IS] and Sastry and Desoer [17] 
are reviewed. 

Finally, in Chapter VII we include, for the sake of completeness, 
some miscellaneous results on the multivariable root locus. These 
include methods for computing graphical bounds on the loci, intersections 
with tne imaginary axis, and other items that might be helpful in plotting 
the multivariable root locus. 

1.2 Notation 

A^(s) represents a scalar polynondal in s. Otherwise, matrices are 
indicated by capital letters, and scalars and vectors are indicated by 
small letters. No underlines are used; whether a quantity is a scalar or 

T 

a vector is clear from the context. A is the transpose of A, DET A the 
determinant of A, TR A the trace of A, Md DIAG [aj...a|^] the n^n mtrix 
with elements ^^...a^ along the main diagonal and zeros elsewhere. ARC z 
is the principal argument of cosqilex variable z, 0(f(x)) the exponent of 



the lowest power of x in a power series expansion of the function f(x), 
and SGN (x] the sign function jf x (one if x is positive, minus one if 
X is negative). 

b^uations, exasqples, figures, lemmas, and theorem: are numbered by 
chapter and position within the chapter; e«iuation (3.17) is the seventeenth 
equation in Chapter III. 

All root locus diagrams are in the complex s plane, imaginary part 
plotted against real part. Open* loop poles are represented by x's, open- 
loop finite zeros by o's, and branch points by triangles. 



CHAPTER II 


Background 


2. 1 Introduction 

In this chapter, we lay the groundwork for the results and analyses 
to follow. First, the basic problem from which the multivariable root 
loci ar^ obtained is described. The characteristic equation is defined 
from tae loop transfer-function matrix, and other equations are defined 
in terns of this equation. These equations are the starting point for 
most of the derivations to follow. The unusual behavior of multivariable 
root loci is accounted for by noting that the root loci are branches of 
an algebraic function. We briefly discuss poles, zeros, branch points, 
break points, and single-point loci, and define these points from the 
characteristic equation. We also define Butterworth patterns and angles 
of root loci. The Newton polygon technique, which gives a series approxi- 
mation to a function of two variables near a zero of the function, is 
described. We will use this technique to obtain results on the angles of 
arrival and departure at zeros and poles. Finally, the SISO (single-input- 
single-output) root locus rules are quickly reviewed for comparison to 
their multivariable generalizations. 

2.2 The Basic Problem 

Consider the feedback configuration shown in Fig. 2.1. and S 2 
are linear multivariable dynamical systems and k is a positive real 
number. It may be shown [2] that the loop transfer-function matrix for 
this configuration is the matrix kG(s), where G(s) is the product of the 
transfer function matrices of $2 and (in that order) . If the feedback 
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loop were broken between and the summation block and signals injected 
at the break, the returned signals at the break would be related to the 
injected signals by minus the loop transfer-function matrix. G(s) is an 
mxm rational matrix function of the complex variable s, and is assianed 
to have full rank and be strictly proper. 

The return-difference matrix for this configuration is I ♦ kG(s), 
and the closed- loop poles are given by the solutions of 

DET [I + kG(s)] s 0. (2.1) 

As k is varied from zero to infinity, the closed-loop poles will vary. The 
plot in the complex plane of the paths swept out by the closed-loop poles 
is the root locus. 

For some results a state-space formulation will be appropriate. The 
system considered is given by 

x = Ax + Bu 

y * Cx 

u = -kKy (2.2) 

where xC r” and y, u C r”*. B and C are assumed to have full rank. Now 
the root locus is the paths swept out by the eigenvalues of the closed-loop 
system matrix 

A^j a A - kBKC (2.3) 

as k varies from zero to infinity. 


I 


IS 


2.3 The Characteristic Equation 
2.3.1 Definitions 

Following the notation of Postlethwaite and MacFarlane [1> 2 , 3], we 
define the characteristic equation 

A(g,s) ^ DET [gl - G(s)l » 0. (2.4) 

From (2.4) we may define two multi-valued functions g(s) and s(g). 
However, these are not ordinary functions of a complex variable, but are 
instead algebraic functions (4). The values of an algebraic function differ 
from those of an ordinary function of a complex variable in that the latter 
form a single analytic function, while the former form a set of analytic 
functions. Each individual function in this set is called a branch of the 
algebraic function. 

The root loci are solutions to 

g(s) ° ~ i ^ positive (2.5) 

and are thus branches of the algebraic function s(g) « 
positive real k [2]. (Single-point loci, to be discussed shortly, are 
omitted.) It is important to note that the multivariable root loci are 
branches of an algebraic function, since this means that their behavior 
can be much more complicated than that of single-input-single-output (SISO) 
root loci. This follows because an algebraic function has as its domain 
a Riemann surface , which consists of several copies, or "sheets," of the 
complex plane that have been "cut" and "stitched together" in such a way 
as to make the function continuous on the surface. A technique for doing 
this is described in MacFarlane and Postlethwaite [1, 2] for A(g,s). 


As k varies, the argument of the algebraic function s(-^) that 
describes the multivariable root loci may pass through a branch point and 
*'jump" from one sheet to another. When this happens, the behavior of the 
root loci may change abruptly *- a locus may swerve, loop around, or act 
in an even more exotic manner (see [5], p. 64). A not uncommon occurrence 
is for a locus on the real axis to abruptly turn around; this behavior is 
discussed in Chapter IV. 

If A(g,s) is reducible (i.e. can be factored over the field of 
rational functions) , there will be several sets of root loci corresponding 
to the several algebraic functions each defined on a different Riemann 
surface. In the extreme case (e.g. G(s) diagonal) where A(g,s) can be 
completely factored into terms of the form (g - g^(s)), gj^(s) a 
rational function of s, the multivariable root locus becomes a super- 
position of m SISO root loci. Reducibility of A(g,s) will not affect any 
of the results of this thesis. 

The following expansion [2] will be used extensively in this thesis: 

4(g,s) - DET (gl - G(s)l = g” - (TR G(s))g®“^ 

♦ iCCprincipal minors of G(s) of order 2)g®“^ - ... 

♦ (-1)® DET G(s) * 0. (2.6) 

Multiplying through by A^^^Cs) , the least common denominator of the nonzero 
principal minors of all orders of G(s) , we obtain 

Hg,s) « Aj„(s)g® + Aj^_j(s)g®"^ + ... ♦ A^(s)g ♦ A^(s) = 0 (2.7) 

where the A^(s) are all polynomials. We can quickly rewrite this as 

'i'(s,g) - B^(g)s‘' ♦ Bj^_j(g)s"’^ ♦ ... ♦ Bg(g) « *(g,s) « 0 (2.8) 
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where the B^(g) are all polynomials and n is the largest degree of the 
A.(s). 

The angle of a root locus s^(k) is the angle a tangent to the locus 
makes with the positive real axis. From Fig. 2.2, 

e = LIM ARG [sjk + 6 k) - s.(k)] 

6 k-K)* ^ 

dSi ds^ 

= LIM^ ARG [-^ 5k] » ARG . (2.9) 

6k-K) 

Note that 

0 » ARG [^] = ARG [^] + ARG (||] 
ds. , ds. 

= ARG [- 5 ^] + ARG [-^] = ARG [- 5 ^]. (2.10) 

A set of loci (s^(k), i = 1 ... r) form an r^^-order Butterworth 

pattern if their asymptotic behavior (as k -^0 or k-^) is such that 

I 

s.(k) 3 (c^k)’' ♦ p^ (2.11) 

where c and o are constants (p_ is called the pivot in the case s-+«), and 
r r r r 

e(Si) = 7 (Qj ♦ 360*i), i * 1 ... r . (2.12) 

Note that this is different from the Butterworth patterns associated with 
optimal root loci, which are the left half-plane portions of larger Butter- 
worth patterns associated with the root square locus. 

2.3.2 Poles, Zeros, Branch and Break Points 
Consider the equation 


*(g»s) * ♦ ... ♦ A^(s) » 0. 


(2.13) 



i 

lj(k+8k) 

8s 

Sj(k) 
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As or we would expect the root loci to either go to infinity or 
approach finite zeros. In fact, setting g * 0 in equation (2.13) gives 
the solutions to 

Aq(s) - 0 (2.14) 

as the points approached by the root loci that stay finite. Similarly, 
making the substitution g • multiplying by (-k)"*, and setting k » 0 
in equation (2.13) gives the solutions to 

\is) » 0 (2.15) 

as the points from which the root loci depart. The question naturally 
arises whether the solutions to equations (2.14) and (2.15) are in fact 
the zeros and poles of G(s). The answer to this question reveals a pro- 
perty unique to the multivariable root locus — the existence of degen- 
erate, single-point loci. 

It is well-known that the pole polynomial p(s) of G(s) is the least 
common denominator of all non-zero minors of all orders of G(s). However, 
Aj^(s) is the led of all non-zero principal minors of all orders of G(s). 
Let e(s) be the led of all non-zero njn-principal minors, with all factors 
common to A^(s) removed. Then 

P(s) - Aj^(s) e(s) . (2.16) 

Since G(s) has full rank, its zero polynomial z(s) is its pole polynomial 
multiplied by the only m^m minor. DET G(s). Then we have 

z(s) » p(s) DET G(s) » e(s) A^jj(s) DET G(s) 

• (-l)\(s) e(s). 


(2.17) 
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Thus the loci that vary with k do in fact depart from the poles of G(s) 
and arrive at the finite teros of G(s) (or at infinity). However, other 
loci start at poles of G(s), remain where they are , and "become" zeros of 
G(s)! It is important to note that these "single-pcint loci ." given by 
solutions to 

e(s) - 0 (2.18) 

are not uncontrollable or unobservable modes of a state-space realization 
of G(s) (decoupling zeros), nor are they poles cancelled when G(s) is formed 
from the transfer function matrices of two separate systems and $2 
(recall Fig. 2.1). They are a phenomenon not found in the SISO case — they 
appear only in the multivariable case. 

However, it is the usual case that e(s) « 1, i.e.. there are no single- 
point loci. And since there is nothing to state about their behavior, they 
will not be considered further. 

Branch points are points where 'l>(g,s) » 0 has a multiple root g^. 

They are associated with unusual behavior of the root locus, since the cuts 
from which a Riemann surface is formed out of several copies of the complex 
plane are made between branch points, or between a branch point and infinity. 
Hence loci can "jump" from one Riemann surface sheet to another, and behave 
strangely, at branch points. 

Break points , short for break-in and break-out points, are points where 
^(g.s) has a multiple root s^. They have the same meaning they do in the 
SISO case. Break points and branch points will be discussed in Chapter IV. 

2.4 The Newton Polygon Technique 

The Newton polygon technique is a graphical device that can be used to 
find a series representation of a function f(x,y) in the vicinity of a 


zero of the function. It plays a central role in this thesis, and under- 
standing it is essential for reading Chapter III. The simple treatment of 
it given here will be sufficient for the purposes of this thesis; for more 
details, see Walker [6]. A theoretical treatment "deriving'* the Newton 
polygon is given first, followed by a step-by-step procedure and an example. 

Consider a function of two variables 


f(x,y) * ♦ Ajy ♦ A^^y" 


( 2 . 19 ) 


with n > 0, Aj^ j* 0 , and A^^ € K{x), where K{x) is the field of all functions 
of X that can be written as a fractional power series in x. 

Puiseux's Theorem [ 4 ] states that K(x} is algebraically closed. This 
means that if 7 is a zero of f(x,y) then 7 can be written as a fractional 
pov T series in x 


CiX 


“1 


C,X 


:2 j.z2 ) 


c^x 




( 2 . 20 ) 


with Cj 0, Zj > 0, Zj^ > 0 for i > 1 (we discard the case 7 « 0, which 
occurs if and only if A^ * 0 ). Puiseux's Theorem is analogous to stating 
that since the field of complex numbers is algebraically closed, any poly- 
nomial with complex coefficients has a complex zero. 

We wish to find possible values of z^ and the values of c^ associated 
with them, i.e. obtain a lowest-order approximation to 7. Substituting 
( 2 . 20 ) in ( 2 . 19 ), we obtain 


^ CjAjX ♦ c^A2X 


, 2z, „ nz, 

^ ♦ cJa x ^ 

1 n 


(Zj+Z2) 2 (2ti+2z2) 


♦ C2A2X 


♦ [CjAjX 

(2z, +Z2) 

♦ 2CjC2A2X ♦ 


« • • 


1 


( 2 . 21 ) 


Since 0, each of the bracketed terms has order strictly greater than 
the order of some unbracketed term. Considering the case of small x, a 
necessary condition for f(x,7) « 0 is that the term of lowest order have 
coefficient zero. Since c^ ^ 0, at least two different terms must have 
the same (lowest) order, and the sum of the coefficients of these terms 
must be zero. Thus there are at least two (and possibly more) indices j 
and k such that 

i jti kz, . iz. 

0(cJAjX n - 0(c*A^x *) i.O(cjA^x ^), i « 0 ... n (2.22) 

where 0(f(x)) is the exponent of the lowest power of x appearing in a 

series expansion of f(x). We may write 

Ai * b^x ♦ (higher order terms) (2.23) 

where * O(A^) (recall A^ CKix)). Then (2.22) becomes 

Sj + jZj * \ <. ♦ iZj, i « 0 ... n (2.24) 

and the sum of the coefficients of the terras of lowest order must be zero: 

^Jb^ . 0 (2. 25) 

summed over the points of the segment giving that value of z^. 

The Newton polygon is a graphical device that yields possible values 
of z^ satisfying (2.24). The polygon for the example to follow is given 
in Fig. 2.3. It is constructed as follows: 

1. Set up a cartesian coordinate system, with u and v axes, and 
plot the n ♦ 1 points P. » (u,v) ■ U,a,), i * 0 ... n, where 
a^ * O(A^) and the A^^ are from (2.19). 
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2. Join to with a convex polygonal arc each of whose vertices 
is a P^ and such that no P^ lies below the arc. This may be done 
quickly by inspection. 

Each segment of the arc defines a line v cu ■ M, where 
is the slope and M is the v-intercept. This line has the property 

V ♦ ou i aj^ ♦ oi, i *0 . . . n. (2.26) 

Let Pj and Pj^ be the endpoints of a segment. Then 

a^ ♦ oj ■ aj^ ♦ ok » V ♦ ou <, a^ ♦ oi, i « 0 ... n. (2.27) 

Clearly Zj ■ o will satisfy (2.24), so that possible values of 
are minus the slopes of the segments of the Newton polygon . 

3. Use (2.25) to find the c^ associated with each z^. For the Zj 

determined by the segment with endpoints P. and P. , we have, for 
the defined in .2.23), J 

c|b^ ♦ » 0. (2.28) 

If there is another point on this segment, we have 

‘^I'^h *^l^k * (2.29) 

This technique may be extended to compute all of the possible z^ and 
c^: see [6]. However, we will not require this extension. 

Example 2.1 We wish to find series approximations to 
f(x»y) ■ (7x“ ♦ 5x ♦ 2)y"* ♦ (8x^ - 6x)y^ ♦ 

♦ (37x^ - llx^)y ♦ (x^° ♦ 4x* - x^) • 0 
in the vicinity of the zero f(0,0) > 0. We identify 

a^ ■ 0, 5^ « 2; a^ > 1, bj ■ -6; Sj ■ 2, b2 • 4 ; 
aj ■ S, bj « -11; “ 6, b^ ■ -1. 

The Newton polygon is shown in Fig. 2.3. Me see that possible values 
for z^ are 

Zj ■ 1, 2 

and the values of c^ associated with then are the solutions to 
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2 

-1 ♦ 4Cj • 0 for 2j * ^ 

2 3 4 

4Cj - 6Cj ♦ 2Cj ■ 0 for ■ I 

which are c. ■ J, -1 for ■ 2 and c. • 1, 2 for z. • 1. So the 
series apprMinations are ^ ^ ^ 

y ■ y ■ y - X. V • 2x. 

2.S The SISO Root Locus 

We quickly review the SISO root locus rules, for comparison with the 
multivariable root locus rules to be given in this thesis. It will be 
seen that some of the SISO rules generalize directly, others less directly, 
and still others not at all. Proofs may be found in any decent elementary 
control theory text (e.g. [7], [8], etc.). 


1. The root locus has n branches, where n is the degree of the 
denominator polynomial of the open- loop transfer function. 

2. All branches of the root locus begin at the open- loop poles. 

m branches terminate at the finite open-loop zeros, where m is 
the degree of the numerator polynomial of the open- loop transfer 
function. The other n - m branches approach infinity along 
asymptotes described in Rule 3. 


3. 


The branches that approach infinity do so along asymptotes 
with angles 


.. i2!Lliy»2:, k . 0, i. : 


(2.30) 


and which intercept the real axis at 


P 




(2.31) 


where the z. are the open- loop zeros and the p. are the open- 
loop poles ^ 


4. A branch of the root locus will lie on the real axis for those 
portions of the real axis that have an odd nusdwr of poles and 
zeros to the right. 


5. The root locus is symmetric with respect to the real axis. 
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The angle of departure from a pole is the sum of the angles of 
all vectors from the zeros to the pole in question, minus the 
sum of the angles of all vectors from the other poles to the 
pole in question, plus 180”. Note that this is 


Vart*i"“= 

+ n * 0, 1, ... k - 1 

th 


(2.32) 


for the angle of departure from a k -order pole p^ of the scalar 
transfer function g(s) . 

The angle of approach to or arrival at a finite zero is the sum 
of the angles of all vectors from the poles to the zero in 
question, minus the sum of the angles of all vectors from the 
other zeros to the zero in question (important : recall the 

definition of "angle" given in Section 2.3.1). Note that this is 


Arrival tCs - 


s = 


"i 


n360* _ _ « , 

+ ' 1| ••• it“ X 

th 


(2.33) 


for the angle of arrival at a k -order zero 


The break-in and break-away points on the real axis may be found 
by solving 

II * 0. (2.34) 

If several branches are approaching and leaving a break point, 
their angles are evenly distributed over 360”. 

If m ^ n - 2, the sum of the closed- loop poles is constant as k 
is varied [8]. 


CHAPTER III 


Angles of Arrival and Departure 


3.1 Introduction 

In this chapter we derive some equations for obtaining the angles of 
arrival (at finite zeros) and the angles of departure (from poles) of 
multivariable root loci. First, the state-space based results of Shaked 
[9] and Thompson [S] are reviewed. Shaked's results employ a spectral 
decomposition; Thompson's results rely on a generalized eigenvalue problem. 
Thus both are computationally arduous. 

Following this, new results are obtained by applying the frequency- 
domain techniques of Postlethwaite [3]. Postlethwaite's approach is to 
obtain a series approximation to the root locus in the vicinity of a 
pole or zero, and then obtain the angles of arrival or departure from this 
series. By investigating how this series is obtained, we derive more 
general results. First, the Smith-NIacMillan form is used to show that 
loci genetically depart from multiple poles, and arrive at multiple zeros, 
in Butterworth patterns whose orders come from the structure indices [12] 
of G(s) at the pole or zero in question. Some equations for the angles of 
arrival and departure are also obtained. 

Next, we use Laurent series expansions of G(s) to derive some simpler 
equations for the angles of arrival and departure. These equations turn out 
to be simple generalizations of the SISO equations. 

Finally, we arrange the coefficients of the Laurent series expansions 
of G(s) mentioned above into Toeplitz matrices, and show how the angles of 
arrival and departure may be obtained from these. We also recall that the 


*• 
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ranks of these Toeplitz matrices are related to the MacMillan orders of 
G(s) at the pole or zero in question, and in this way all the results of 
this chapter are tied together. 


3.2 Review of State-Space Results 

The results described in this section were first derived by Shaked [9] 
and Thompson [5] . They are included here for comparison to the frequency 
domain results of the following sections. Also, when a system is speci- 
fied by its (A,B,C) matrices these equations may be simpler, since compu- 
tation of G(s) * C(sl - A)“^B may be very difficult. 


Theorem 3.1 The angles of departure are given by 


0. ^ »ARG [-v!bkCu.] 

depart ^ 


(3.1) 


T 

where Uj^ and Vj^ are the right and left eigenvectors associated with 
the open-loop pole (an eigenvalue of A) considered, and Uj^ and v} 
have been normalized so that vTuj^ * 1 . 

The angles of arrival are given by 


6 . = ARG (wTk‘^x.1 

arrival '• i i^ 

T T 

where w^ and solve, for some y^ and q^. 


A - z^l 

-C 



f "" 

•H 

X 

1 


— 1 

•H 

X 

1 


= 0 


(3.2) 


(3.3) 



(3.4) 


T 

where is the finite zero considered, and Xj^ and have been nor- 
malized so that wTx^ « 1. 


Following Shaked [9], one may prove (3.1) by recalling that the closed- 
loop system matrix is 


A ,» A - kBKC 
cl 


(3.S) 


and by letting its spectral decomposition be 


DIAG [s.l [Uj ... u^l (3.6) 

where the are the closed- loop eigenvalues with eigenvectors and left 
T 

eigenvectors v^^, all of which are functions of k. We have 



''ci“i • n“i 


C3-7) 


T 

and differentiating with respect to k and multiplying by v. , one gets 


,t"c1 T, 

'^i dk ‘^cl dk 


i T''"i 
Tr- + s.v. -TT- 
dk 1 X dk 


(3.8) 


so that 


ds. 

X 

■ar 


T ^cl 


vT BKCu.. 

X X 


(3.9) 


Letting k 3 0 and using equation (2.9), (3.1) is proved. Thonq>son [5] 
derives this same result from a generalized eigenvalue problem, and his 
method is preferable for proving (3.2). Indeed, in [5] Thompson points 
out several errors in Shaked's paper. 

T 

It should be evident that s^, u^^, and v^ solve the generalized 
eigenvalue problem 




0 




(3.10) 
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[■: ■:] 


A - s.I 


B 




(3.11) 


T 

for some and x^. 

T T 

Differentiating (3.10) with respect to k and multiplying by [v^ w^] 
gives> using (3.11), 



(3.12) 


which can be rearranged into 




T 


w. 

1 



X. . 
1 


(3.13) 


Letting using equation (2.9), and noting that ARC [-—r] = 0 proves 

k^ 

(3.2) > Note that Rosenbrock's [10] definition of finite zeros guarantees 
a non-zero solution to (3.3) and (3.4). This concludes the proof of 
Theorem 3.1. 


3.3 Results from the Sroith-MacMillan Form 

We use the Smith-MacMillan form of G(s) to obtain equations for the 
angles of departure and arrival using the polynomials A^(s) of the charac- 
teristic equation, and to characterize the loci departing from and arriving 
at multiple poles and zeros. For simplicity and clarity we consider first 
the case of first-order poles and zeros, exhibiting the methodology, and 
‘then proceed to the far more complex case of multiple poles and zeros. In 
both cases we consider angles of departure first, and then the angles of 


arrival. 


3.3.1 First-Order Poles and Zeros 
We start with 

«(g,s) « Aj^(s)g® + ♦ ... + A^(s) • 0. (5.14) 

Since we are first interested in angles of departure, i.e. the behavior of 
the root loci as k-*”0 and g-*--«, we make the substitution g * and 


multiply by (-k)®, yielding 

AoCs)(-l)®k® ♦ AjCs)(-l)®‘^k®"^ ♦ ... - Ajjj_j(s)k (3.1S) 

+ A (s) * 0 . 
m 

Let be a first-order pole, and define the following: 

s = s - p^ (3.16) 

A^(s) » Ajjj(s + p^) a ajS + a2S^ + ... (3.17) 

Pi^ ' *^0 **1® * * •*' (3.18) 


where the right sides of (3.17) and (3.18) are finite polynomials. Note 

that the constant term of (3.17) is zero, since A (p.) « 0. 

m i 

By applying the Newton polygon technique to (3.17) and (3.18), we may 
obtain a series approximation to the locus in the vicinity of the pole p^, 
and obtain the angle of departure at once. In particular, if b^ and a^ are 
both non-zero, the Newton polygon will be as in Fig. 3.1. We then have 

k a c5 as (S,k)-*>(0,0) (3.19) 


where c solves 



a, ♦ cb « 0. 
1 o 


(3.20) 


This gives 


s s --- k as (s,k)-*> (0,0) 


and, using the definition of angle, 

b. 


depart 


= ARC 


(3.21) 


(3.22) 


We will refer to the case where a, and b are both non-zero as the 
generic case. The word "generic" is used to describe a property of a 
finite set of parameters which holds for all values of these parameters 
except those satisfying a finite number of polynomial equations [21]. Thus 
a generic property "almost always holds." The property that a^ and b^ are 
both non- zero is indeed generic, since it will be shown shortly that a^ 
is always non-zero if p^ is a first-order pole, and that b^ is zero only if 
p^ is a root of both Ajjj(s) and Ajj|_^(s). This can occur only if the coeffi- 
cients of A^(s) and A^jj_j(s) satisfy a polynomial equation obtained by 
setting the resultant (see Appendix) of A^(s) and A^_^(s) equal to zero. 

Now write G(s) in the Smith-MacMillan form 

rni(s) n (s)-i 

G(s) «U(s) DIAGj^^^, ... , Jil^JV(s) (3.23) 

where U(s) and V(s) are unimodular (have constant determinants) and where 
nj^(s) ! n^(s) and d^(s) I dj^(s) for 1 i i i j i m. (Recall that G(s) is 
mxm and has full rank.) Also, let 


V(s)U(s) » W(s) « [w.j(s)l . 


(3.24) 


We may write 
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A(g,s) - DET [gl - G(s)l » DET [gl - U(s) DIAG 
DET [gl - DIAG|:ji77T|V(s)U(s)] 




*n^(s)* 


V(s)] 

(3.2S) 


and we may obtain <l>(g,s) from the right side of (3.25). 

We now state and prove the following theorem: 

Theorem 3.2 If and only if t 0, then we have 


depart * 


[-^I.ARG 

\-l"’ I 

<1 , X * ® “ P: 

L ij 



(3.26) 


(Note that this corresponds to the "generic case" discussed earlier.) 

We prove Theorem 3.2 by showing that a^ and b^ are non-zero if and 
only if Wjj(p^) is non-zero. Then the theorem follows immediately from 
(3.22), while the expressions for a^ and b^ follow immediately from (3.17) 
and (3.18). 

We observe first that, since p^ is a single pole, 

Ajj(s) » (s - P^)\(s) * “ *1* + ••• (3*27) 

where \(pj^) i 0. Dividing by s and setting s » p^^ immediately gives 




(3.28) 


As for b^, note that 
0 


and since p^ is a single pole, we may write 


(3.29) 


d^(8) ■ (s - Pj^)dj(s) 


(3.30a) 


di(P^)i d2(p^)f ... djij(P^) ^ 0 


(3.30b) 


so that (3.29) becomes 


n,(s) 


TTT “u‘” * V‘)2 


4.0) 


Hj^Cs) 

3~7aT 1 (s) * (3.31) 

i-2 1'“*^ 


Since n^(s) and d^(s) are relatively prime, n^^(p^) is non- zero. Therefore, 
by combining (3.18), (3.30b), and (3.31), we find that 

”o ■ -Vl'Pi> • 


4, CP,) 


» n.(p.) 

* VPPE 3jTpIT"jj(Pi)- 


(3.32) 


Since the second term vanishes, b^ is non-zero if and onl)' if w^CP^^) is 
non- zero. This proves Theorem 3.2. 

Let us now consider the angle of arrival at a finite zero z^. Since 
we are interested in the behavior of loci as or g-^0, we may work with 

«(gpS) » ♦ ••• ♦ Aj(s)g ♦ A^(s) • 0. (3.33) 

We now make new definitions 

5 » s - Zj^ (3.34) 

A^(s) « A^(S + Zj^) » a^i ♦ a 2 i^ ♦ ... (5.35) 

A^(s) » Aj(5 ♦ z^) » bg ♦ b^i ♦ b 2 i^ ♦ 




(3.36) 


As before, we consider the generic case when and are non-zero, for 
which the Newton polygon is illustrated in Fig. 3.2. The resulting cal- 
culations parallel (3.19)-(3.22) (using g instead of k) , and the following 
theorem should not be surprising: 


Theorem 3.3 The angle of arrival at a first-order zero z, is given 

iy ^ 


e 


arrival 


. argI-^I.arg I*. 


(3.37a) 


if and only if 

Wii(Zi) 


DET 


•• **1 m-l^*i^ 

• • 

• • 

• • 

Vl ••• %-l 


* 0 . 


(3.37b) 


However, the proof of Theorem 3.3 is more difficult than the proof of 
Theorem 3.2. It requires the following lemma, which will be used exten- 
sively in the next section. 


Lemma 3.1 (Binet-Cauchy Theorem) Define the following notation for 
pxp minors of an nxn matrix A: 


••• SY 

v> ■■■ V 


DET 


*i i ^i i 


L\J 


a. 


pn 


Vpj 


, lii, < i, < ... < i in 

12 p 

lijj < < ... < jpin 


Then 




summed over all possible l<.k^<k 2< ... <k^in. 

This standard result is given and proved in Gantmacher [11]. 
To show that a^ 0, write 


(3.38) 



(3.39) 


\(s) ■ ys)(-l)"'(DET G(s)J 

- Ajjj(s)(-1)®(det u(s)Hdet V(s)]^n^ 

and note that since is a single zero we may write 

nj||(s) • (s - 2i)n„(*) (3.40a) 

nj(z^)« ... ^ (3.40b) 

Also note that, excluding single-point loci, 

m 

V*) • " (3.41) 

i*l 

and recall that U(s) and V(s) are unimodular, so their determinants are 
non-zera constants. Then we have 

Aq(s) • (-D^tDET U(s)HDET V(s)H^n^ n^(s)]h^(s)i (3.42) 

and by comparing this to (3.35), dividing by i, setting s * z^, and using 
(3.40b), we see that a^ is non-zero. As for b^, using Lemma 3.1 on (3.25) 
gives 


Aj(s) . A„(s )(-l)®'^ £ (principal minors order m-1 of G(s)J 


,«/ "<(»)\ / corresponding \ 

• \,(s)(-l)"‘ principal 1 

Vr**^m-1 J / of W(s) / (3.43) 

since all the non-principal minors of the Soith-MacMillan form are zero. 


However, since ■ 0, 


*>0 ■ 

■ A. 


„ , m-1 n,( 2 .) f ••• '*1 


44 ) 


and recalling (3.40b) and (3.41) shoMS that is non-zero if and only if 
(3.37b) holds. This proves Theorem 3.3. 

In Section 3.4, we will obtain some alternative equations for the 
angles directly in terms of G(s) and G~^(s). 


3.3.2 Multiple Poles and Zeros 

Before we investigate the angles of arrival and departure for multiple 
poles and zeros, some discussion will be necessary on exactly what is meant 
by a multiple pole or zero. This is not a trivial matter; in the multi- 
variable case, it is possible for G(s) to have a pole and a zero at the 
same location, or several poles and zeros of various orders all at the same 
location. We now make several definitions that will clarify matters and 
make the analyses to follow as straightforward as possible. 


.th 


Definition The pole p, is a k -order pole if the exponent of 
(s • p^) in the pole polynomial of G(s) is k. 


n^(s) 


Definition Let the 

and let p. be a k^^-order pole. Let k. be the largest integer such 
that ^ ^ 


Smith-MacMillan form of G(s) be DIAG [^~^1 • 


k. 


(s - p^) J ! dj(s), j ■ 1 ... m. 


(3.45) 


Then the ik. ) are the structure indices [12] associated with the pole 


'i- 




Analogous definitions are made for k^-order zeros and structure 
indices of zeros. Note that p^ may be a k^^-order pole and also a zero. 
Also note that p^ may have one set of structure indices as a pole and 
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another set of structure indices as a zero. 

k. a k 
J 

(2) For a polQ, k > k- i k, i ... i k 

i j in 

For zero, the ordering is opposite 

(3) For a first-order pole, the structure indices are 

[1, 0, ... 0j« 

Example 3/ ; Suppose that G(s) has as its Smith-MacMillan form 

1 , s - 1 ^ cs.-.ir] ^ 

(s + 2)^(s ♦ 1)^* (s ♦ 2)(s + 1)^’ s + 1 J 

Then the structure indices are 

Pole at -2: [3, 1, 0] 

Pole at -1: [2, 2, 1] 

Zero at 1: [0,1,2]. 

It should be noted that this definition of ’'structure indices" is not 
the same as the definition given by Verghese and Kailath in [12], although 
the definitions are closely related. The difference may be illustrated as 
follows. Let q be both a pole and a zero of G(s), with Verghese-Kailath 
structure indices [a,, a_, ... a., 0 ... 0, -o., ... -a ]. What we shall 
do here is separate the polar structure of q from its zero structure. From 
this point of view, the structure indices of G(s) at the pole of q are 
given by [Oj, a^, ... a^, 0 ... 0] and the striicture indices of G(s) at the 

zero q are given by [0 ... 0, a., ... a ]. 

J m 

The motivation for doing this is that the polar nature of q has no 
effect on the angles of arrival at q, and the zero nature of q has no 
effect on the angles of departure at q. So we may consider q as consisting 
of a pole and a zero which just happen to be at the same location, but may 
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E 


be treated separately. 


Having defined terms and notations, we proceed now to investigate the 
angles of arrival and departure for multiple poles and zeros. The metho- 
dology and results employed in this section will be the same as those of 
the last section, to which extensive reference will be made. Once again, 
the angles of departure are treated first, and the angles of arrival 
later. 

Let p^ be a k^^-order pole with structure indices [kj, k 2 , ... k^j^]. 

We require series approximations for all of the loci departing from p^, 
and thus we must determine the form of the Newton polygon from the structure 
indices. As before, define 


S = s - p^ (3.46) 

2 

A^(s) = Ajjj(5 + p^) - * ajS + a^s + ... (3.47) 

“ ’Vl^^ * ?i^ * ‘’o * ‘’l^ ^ 

A^_2(s) * ~ * *^1* ‘^2^” * (3.49) 

A^(s) . (s - Pi)\.(s), V?i) # 0. (3.50) 


From (3.50) and the same argument as used on (3.27), 

a^ — a^^ — ... — a^^ ^ ® 0, a^^ ^0. (3.51) 

We define 

k . ^ 

dj(s) » (s - p^) ^ dj(s), dj(p^) 0, j * 1 ... m (3.52) 

and recall (3.29) to write (corresponding to (3.31)) 


(3. S3) 


(k-k.). n.(s) 

• Pi’ ''m'®’ "ll‘” 


ra (k-kj, n.csj 

*E (S-Pj) ’A^(s)^W Cs). 
J*2 dj(sj 


n.(s) 


Comparing (3. S3) to (3.48), dividing by increasing powers of 5 and setting 
s = p^ (a familiar procedure by now), we get 


. bj . ... -b^.k^.i - 0. b,.^^ t 0 


(3.S4) 


where we have assumed (1) Cp^) ^ 0» (2) kj k 2 » If kj » k 2 » we require 
instead the assumption 


0 (5.SS) 


/ni(Pi) n 2 (p.) ,\ 

\ k * ^ "'22^Pi^ J 

>'-''1 ^ ‘ Vdi(Pi) ‘‘zCPi) ' 

which we would still expect to be true in general. 

Now we must consider k -,(s). Again using Lemma 3.1, we may gener- 

m-A 

alize (3.43) to 


A (s) * Ajjj(s) (-1)™’^ 2 [principal minors order m-p of G(s)] 



n 

••• Vp 



corresponding 
principal 
minor of W(s) 


) 


» 


psi ...m. 


(3.S6) 


Taking p = m - 2 and some reflection leads to 


V2^"i • 


(k-k.-k,), n, (s)n-(s) 

(s - Pi) ^ ^ A^(s) . DET 




Wii(Pi) Wj2(Pi) 

dj(s)d2(s) ["’2ltPi^ '‘'22^Pi^, 

/ corresponding 2x2 A 
sj/ \principal minor of W(s)J (3.S7) 


V V 
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(recall kj k 2 i i i k^j^). Since n^Cs) and dj(s) are relatively 
prime, and is zero, n^Cp^^) is non-zero. A similar argument applies 

for n 2 (P|) if k, is non-zero. 

Comparing (3.57) to (3.49) and once again dividing by increasing 
powers of s and setting s ” P^^i we get 


c^ a c, « ... 
o 1 


Vkj-k2-l * ° 


(3.58) 


where we have assumed k. / k. and 

2 0 

Wii(Pi) Wj2(?i) 
W2i(Pi) «22(Pi) 


DET 


* 0. 


(3.59) 


If k 2 a k. we must satisfy instead an equation analogous to (3.55). 

We may repeat this argument for p a ni-3, m-4, ... , and obtain more 
equations analogous to (3.51), (3.54), and (3.5S), with analogous assump- 
tions. The Newton polygon will take the form given in Fig. 3.5, and we 
can now prove the following theorem: 


Theorem 3.4 The loci departing from a k^^-order pole depart gener- 
ically in Butterworth patterns whose orders are the non-zero structure 
indices of G(s) at the pole. For a pole p^ with non-zero structure 

indices [k^, k 2 , ... k^] the angles of departure are: 


r (k-k.) 
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’depart, 1 k^ 


ARC 
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CFkp'Vi'®’ 


,(k) 
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ds(k) m 


A„(s) 


Pi 




n360* 



n*0, 1, ... 
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(3.60a) 




3 * y * 

depart ,2 k. 


ARC 


(k-kj-k2) 


ds 


(k-kj-k2) 




(k-kj) 


S = p. 


L ds 


(k-kj) 


♦ il^£_ ^ n=0, 1, ... k2 - 1 


3 * T*” 

depart , r k 


ARC 




(lO 


Pi 


ds 


(kp Vr+1^^^ 


^ n360“ , 

t n ® Oj If ... 


if the following assumptions are met; 

(1) kj / k2 >^ ... / kj. 

w^j(Pi) ... Wjj(p.) 

• • 

• • 

• • 

Kjl(Pi) ... 


(2) DET 


^0| J “ 1 ... r. 


We prove Theorem 3.4 by applying the Newton polygon technique, 
the are ordered, the polygonal arc drawn in Fig. 3.3 is convex, 
indeed the Newton polygon. If the k^ are unequal, we may write 

k s c^S , i » 1 ... r, as (5,k) -*■ (0,0) 


where the c. solve 
1 




(3.60b) 

(3.60c) 

(3.61a) 

(3.61b) 

Since 
and is 

(3.62) 


(3.63) 


where is defined as 


(k-kj-.. .-k.) 

Ajjj^i(s) a 3^s +... ,i*l...r (3.64) 

(note 3^ » 3^ ■ , and 3j ■ )• From (3.62) and (3.63), we 

get 


1/k, 


x.k 

1 


(3.65) 


where solves 




i 


1 ... r . 


(3.66) 


The k^ solutions to (3.66) are equally spaced in the complex plane by 
angles of (360 /k^)**; hence their angles are those of a Butterworth pattern. 
This fact with (3.65) proves that the departing loci form Butterworth 
patterns. Applying the definition of angle of a locus to (3.65) and noting 


that 


®i 


1 

(k - kj^ - ... - kj^) ! 


(k - k. - ... - kj 
d ^ ^ 

(k - k, - ... - k.) Vi^*^ 
ds ^ ^ 


‘s » p^ 
(3.67) 


yields (3.60a, b, c) and concludes the proof of Theorem 3.4. 

Note that if any two non-zero k^ are equal then three points in Fig. 
3.3 will be collinear, and (3.63) will be replaced by an equation o£ the 
form 


®i-l ^ ""i^i * ‘^r^^i+i 


(3.68) 


which is now quadratic in c^. In this case, we get two Butterworth patterns 
of order k^ with different principal angles, whose computation from the 


A^(s) is now much more difficult. 


As for the angles of arrival at a multiple finite zero, we may use a 
procedure analogous to that used for the angles of departure to obtain an 
analogous result. The groundwork was laid in the last section for first* 
order zeros; the only difference is that we now use p « 1, 2, ... r in 
(3.56) rather than p > m - 1, m * 2, ... m - r. The result is 


Theorem 3.5 The loci arriving at a -order zero arrive generically 
in Lutterworth patterns whose orders are the non- zero structure 
indices of G(s) at the zero. For a zero z^ with non- zero structure 

indices [k^, k_ ,, ... k„ the angles of arrival are: 
in ni*i lu-r 


arrival , 1 “ Tc 

m 


(k-kjjj) ^1^®^ 




'S » 2 


««• k " i 
m 


1 ^r+l^*^ 

^arrival, r + 1 * I — (k ) 

m-r m-r-' 

d 


rk — T 


S * 2 . 
1 


^ n360” , ft 1 V 1 

* Ic ~ »n»0, 1, 

m-r 


if the following assumptions are met: 

^m ^ *Sn-l ^ ^ Vr 


(3.70) 


«ll(Pi) «ij(Pi) 

(2) DET : : j* 0, j «m-l,m-2. . 

'^ji’cPi) ... 


. . m - r - 1 . 


(3.71) 


(For notational convenience, we let r be one less than the number of non- 
zero indices.) 

It is unfortunate that the assumptions (3. alb) and (3.71) are so 
difficult to verify, depending as they do on the unimodular matrices 
U(s) and V(s) that bring G(s) to the Smith-MacMillan form. Except for one 
note in the next section, it has not been possible to interpret these 
conditions or come up with easier ways to determine whether or not they 
are fulfilled. In any case, it is felt thnt the results of this section 
cast considerable light on the angles of root loci departing from or 
arriving at multiple poles or zeros. 

We end this section with an illustrative example: 


Example 3.2 We wish to find the rsigles of arrival and departure for 
the root locus of 


G(s) 


1 

s'* + 4s^ ♦ 8s^ + 8s ♦ 4 

I s^ + 8s ♦ 17 
s^ + 9s^ + 2Ss f 17 



s^ lOs^ ♦ 33s ♦ 34 
♦ 21s^ ♦ 78s" ♦ 117s ♦ 68 


We have 


^(g.s) » g^ - TR G(s) ♦ DET G(s) » 0 

and it is straightforward to compute h(g,s) and multiply through by 
the least common denominator to get 

♦(g.s) • (s^ ♦ 6s^ ♦ 18s^ ♦ 32 s^ ♦ 36s^ ♦ 24s ♦ 8)g^ 

- (2s^ ♦ 25s* ♦ 125s^ ♦ 32Ss^ ♦ 493s^ ♦ 420s ♦ 170)g 

♦ (s^ ♦ 16*i^ ♦ 98s^ ♦ 272s ♦ 289) « 0 . 

The open- loop poles are found by solving 

A2 (s) * s^ ♦ 6s^ ♦ Ids'* ♦ 32s^ ♦ 36s^ ♦ 24s ♦ 8 • 0 


which yields the third-order poles -1 ♦ j and -I - j. The finite 
zeros are found by solving 

A (s) » s^ ♦ 16s^ ♦ 98s^ * 272s ♦ 289 -0 
o 


which yields the second^order zeros -4 ♦ j and -4 - j. The structure 
Indices may be found by constructing the Smith-MacMillan form of G(s) 
(another method will be described later): 



0 

s^ ♦ 8s ♦ 17 
s^ ♦ 2s ♦ 2 . 


from which we see that the structure indices for the poles -1 i j are 
[2, 1), and those for the zeros -4 ± j are [l, 1]. Since we have 
U(s) and V(s) , we may quickly confirm that the W(s) assumptions are 
upheld. 

Thus loci depart from the pole •! ♦ j in some second-order and 
first-order Butterworth patterns with angles 


9 3 4> 

depart , 1 2 


ARG 






Lds 


(3) "2 


A,(s) 


8 » p. 


♦ nl80®, n » 0, 1 


1 r 12s^ ♦ I25s^ ♦ SOOs^ ♦ 9758^ ♦ 986s ♦ 420 i 

_ Arg I s 1 s » -1 ♦ j 

^ L 120s^ + 360s^ ♦ 432s +192 » * jj 


61.8", 241.8* 


and angle 


depart , 2 


ARG 


*,t») 


« ARG^ 

» 33.7*. 


s'* ♦ 16S^ ♦ 98s^ ♦ 

. 125^ ♦ 12Ss^ ♦ 5005^ + 9 


s • Pi 

« 

272s ♦ 289 


97Ss^ ♦ 986s ♦ 420 


. -1 ♦ j j 


By symmetry, the angles of departure from the pole -1 • j will be 
-61.8*, 118.2*, and -33.7*. 

We know that the loci will arrive at the finite zero -4 j in 
two first-order Butterworth patterns, but since the structure 
indices are equal, the computation of the angles of arrival would be 
much more difficult, and is not attempted here. 

3.4 Results from Laurent Series 

We now use a different methodology to come up with much simpler 
equations for the case of first-order and certain higher-order poles and 
zeros. These turn out to be nice generalizations of the SISO equations. 
Let p^ be a first-order pole of G(s) , and let the Laurent expansion 


of G(s) at Pj^ be 

G(s) « g l^ ' p7 ^-1 ^0 ^ " ^i^^l * “ Pi^^^2 * *** 

Then taking the trace of both sides, we have 

TR G(s) « ^ - j - p - TR G_J ♦ TR Gq ♦ (s - TR + ... (3.73) 

Recall (3.16-18): 

5 • s - p. (3.74) 

A^(s) ■ \(5 + pj^) * aj^S ♦ aji^ ♦ ... (3.75) 

-Aji-i(s) « ♦ Pj^) ■ * ^1* ■*’ * *** (3.76) 

Now multiply (3.75) and (3.73): 

» (a^ TR G_j) ♦ i(aj TR G^ ♦ aj TR G,^) ♦ ... (3.77) 


Equate coefficients of (3.76) and (3.77): 


bo . aj TRG.^. 


(3.78) 


We know from (3.28) that is non-tero, and recalling (3.22), 

« d.p«t • *«<= t~i 
we have proved Theorem 3.6: 


(3.79) 


Theorem 3.6 Let the Laurent expansion of G(s) at a first-order pole 
Pj^ be (3.72). Then, if TR G^j i 0, the angle of departure from pj^ is 


® depart * ^RG 


t-TR G^j] - ARG t-(8 - pp TR G(s)|j,p^]. (3.80) 


Simple as this result is, it is a striking generalization of the SISO 
root locus equation (2.32) for computing the angle of departure. The only 
difference is that in (3.80) the trace of the transfer function matrix is 
used, whereas in (2.32) the (scalar) transfer function itself is used. 

The condition TR G_^ #0 will hold in general; in fact, it is easy to 
show that TR G_j f Q if and only if w^j(pj^) i 0, which is the condition 
that is needed in Theorem 3.2 in order to use (3.26) to compute the angle 
of departure from the A^(s). Recall the Smith-MacMillan form of G(s) 


G(s) « U(s) DIAG 




and use (3.30b) (properties of the d^(s)) to show that 


(3.81) 


G,i • (s - Pi)G(s)lj . 

[ ni(Pi) 
di(Pi) 


, 0 ... 0 V(p.). 


(3.82) 


Then using the commutative property of the trace, one has 


TR G 


TR 



"l(Pt) 

djCPi) 


ni(Pj) 


Wll(Pi) 


0 



(3.83) 


and the result follows. 

In the unusual case where TR is zero, something rather interesting 
happens. Ne now state and prove 


Theorem 3.7 Let the Laurent expansion of G(s) at a 
p^ be (3.72). Then if 

first-order pole 

(1) TR G_j • 0 

(3.84) 

(2) TR [G,j G^] J* 0 

(3.85) 

the angle of departure from p. is 

f* ^ T 


»d.p«t • "><= [™ "=-l =olJ 

(3.86) 

and the locus departs as k" (a 1-order departure) . 



If TR G_^ is zero, is zero, and the Newton polygon for this situa- 
tion, given in Fig. 3.4, shows that we must consider 

V:‘*) • *.(•) S 

Recalling the Laurent expansion of G(s) , we see that 

L /principal minors\ 1 ^ /principal minors\ 

yorder 2 of G(s) J _ p ^2 " ^order 2 of G^j J 

1 « /principal minors with one column 

s - p^ ^ yfrom G^^ and one column from G^ 

and the first term must be zero, or p^ would be a second-order pole. 

Denote the elements of G , and G. as 

-1 o 


... (3.88) 



54 






(3.89) 


and observe that 



principal minors with one column 

from G , and one column from 

-1 o 



0 -1 0 \ 
" ®ji ®ij ) 




» (TR G_p(TR G^) - TR [G_^ G^] = -TR [G_^ G^j. 


(3.90) 


Multiplying (3.88) by the polynomial for A^^(s) (3.75) and using (3.87) 
and (3.90), we get 


TR [G_j Gq] -k s(-a2 TR [G_^ G^] + . . . ) (3.91) 

so that if TR [G_^ G^] is non-zero, the Newton polygon will be as in Fig. 
3.4. We then have 

k s c§^ (3.92) 


where c solves 

a^ - a^c^ TR [G_j^ G^] =0. (3.93) 

Since a^^ is known to be non-zero from (3.28), we have 

S s TR [G_^ G^]k^ (3.94) 

and the theorem follows. 

Unfortunately, the Laurent series methodology does not lend itself 
well to the case of multiple poles. In general, the results from the 
Smith-MacMillan form must be used. However, for a certain class of 


multiple poles, the following result applies: 


fh 

Definition A k -order pole is said to be simple if its structure 
indices are [k, 0 ... 0]. 

Theorem 3.8 Let the Laurent expansion of G(s) at a k^*'-order pole 
Pi be 


G(s) 


(s “ Pi) 


^ ■ ' . G ^ ■ G I ♦ G ... 

K f ^ \ O 


(3.95) 


(s - Pi) 


Then if TR G_^ ^ 0, the pole is simple and the angles of departure 
from Pi are 


1 . 


n360‘ 


depart " ^ " • 0, 1, ... k - 1 


^ ARC (-Cs - Pj)'' TR 0(8) Ij , p i 
nSoO® 


T 


*, n ®0, 1, ... k — 1* 


(3.96) 


The proof of Theorem 5.8 follows that of Theorem 5.6. Recalling 
(5.51) for a k^^-order pole we have that 


a » a, = ...» a, , » 0, a. / 0. 

0 1 K-l K 


(3.97) 


Taking the trace of (3.95) and multiplying by \(s), we get 

a (aj, TR G_J^) ♦ S(aj^^j TR G_j^ ♦ aj^ TR G_j.^i) ♦ ... 


(3.98) 


Equating coefficients with (3.76) gives 


b a a. TR G 1 . 
ok -k 


(3.99) 


Now, if TR G j. is non-:ero, b^ is non-:ero, and the Newton polygon will 
look like Fig. 3.5. Comparing Fig. 3.5 with Fig. 5.3 shows that Pi must be 



Figure 5.5 

)Iygoti for the Generic Case 
Simple Higher-Order Pole 


simple. We have 
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* depart 

where x solves 

. b 

» - - 2 . a _TR G . (3.101) 

h. 

as before (see (3.64-66)). (3.100) and the k solutions to (3.101) yield 

(3.96) > and the theorem is proved. 

If TR is zero, p^ is not a simple pole and we must go back to the 
Smith-MacMillan form results. The reason that the SISO angle of departure 
equation (2.32) does not always generalize to (3.96) is that in the SISO 
case all higher-order poles are necessarily simple, while in the multi- 
variable case only some are. 

The Laurent series methodology might at first seem inapplicable to 
the computation of angles of arrival, since these depend on the sum of 
principal minors of order m - 1 of G(s), which, unlike the trace, is not a 
linear function. However, recall that 


G'^(s) 


ADJ G(s) 

DEY-^s7 


(3.102) 


where ADJ G(s) is the transpose of the matrix of cofactors of G(s). So the 
main diagonal elements of G~^(s) are principal minors of order m - 1 of 
G(s) divided by DET G(s), and we have 


A^(s) 


. . .^/principal minors of\. .xm-l 
m'-*^^yorder m - 1 of G(s) ^ 


* Ajjj(s) DET G(s) 

« -Aq(s) TR G"^(s). 


^mai 


main diagonal elements 


of G'^(s))(- 


1 ) 


m-1 


I 


i 


i ^ 



(3.103) 
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We once again define, for a k^^-order finite zero 


S * S 


(3.104) 




Ao(s) 

• A^(5 2i) » a^ ♦ a, S + a-S“ + ... 

0 10 1 2 

(3.105) 

AjCs) 

1 

a A (S ♦ z.) a b^ ♦ b.S ♦ b,S" ♦ ... 
1 1 0 1 2 

(3.106) 


where we know (see (3.39-42)) 

a ... a * 0» 0 * (3.107) 

Since is a k^^'-order zero of G(s), it is a k^^-order pole of G"^(s), 
and we may write a Laurent expansion 

G"^(s) r H .■*•... + 1 M , + H + ... (3.108) 

(s - z^)‘^ ■ (s - r^) 

Taking the trace of (5.108), multiplying by (3.105), using (3.107), and 
equating coefficients with (3.106), we get 


bo a -a^ TR H., 


(3.109) 


and following (3.99-101), we have proved 


Theorem 3.9 Let z^ be a k^^-order zero of G(s). Then if TR ^ 0, 
where is defined by (3.108), the angles of arrival at z^ are 


0 


arrival 


^ ARG [(s - z.)*" TR G’^(s) I s = 2. ] + 
n a 0, 1, ... k - 1. (3.110) 


As with the angles of departure, TR being non- zero corresponds to z^ 
being a simple zero. If z^ is non-simple, the Smith-MacMillan form results 


must be used 
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We can simplify (3.110) for the case of 2 x 2 G(s), by noting that in 
this case we have 




(3.111) 


so that (3.110) becomes 


arrival “ Ic 


[ 


(s - TR G(s) 
DET G(s) 


] n360* 


n»0, 1, ...k-l« 


(3.112) 


We conclude this section with a simple example: 


Example 3.3 [3] We wish to find the angles of arrival and departure 
for 


G(s) = 


s'* + 5s^ - 2s^ - 44s + 40 


3s^ ♦ 4s" - lS6s + 464 
s^ ♦ 79s^ + 44s - 868 


8s - 24s ♦ 16 
-4s^ - 4s^ + 40s - 52 


From the characteristic equation (or other means) , it is ascertained 
that G(s) has first-order poles at 1, 2, -4 + 2j, and -4 - 2 j , and 
first-order finite zeros at 1 + j and 1 - j. Using Theorem 3.6, we 
have 


Pole at s s 1 


« depart = ^RG [-(s - 1) TR G(s) ^ , j] 

• or r s^ ♦ 116s - 432 , 

* [(S - 2)(s + 4 + 2j)(s ♦ 4 - 2j) > s = 1 

Pole at s « 2 

^ f s^ ♦ 116s - 432 , 

° depart “ L(s - l)(s + 4 + 2j)(s + 4 - 2j) • s » 2 

= 180 * 
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Pole at s ■ -4 ■» 2j 

A Aor r ♦ 116s - 432 , 

depart * ^(s - l)(s - 2) (s ■*•4 + 2j) ' s * -4 ♦ 2j 

» 110.9* 


Pole at s ■ -4 - 2j 

® depart “ -110-9“ synunetry . 
Using Theorem 3.9 and (3.112), we have 


Zero at s » 1 j 



Zero at s « 1 - j 

* -69.6“ by symmetry. 


These results agree with those of Postlethwaite [3], from which this 
example is taken. 

3.5 Results from Toeplitz Matrices 

We now take the Laurent series coefficients from the last section and 
arrange them into Toeplitz matrices from which the angles of arrival and 
departure can be obtained. We also relate these results to the preceding 
results, showing how all of these results relate to each other. The 
results that follow are due to Sastry and Desoer [17] and Levy [19]. 


Definition A matrix A has simple null structure if, in the Jordan 
form o^ the zero eigenvalues are all contained in Jordan blocks 
of order one. 


Note that if a matrix has simple null structure, then its rank 


is equal to the number of non-zero eigenvalues. 


Theorem 3.10 Let the Laurent series expansion of G(s) at an n^^-order 
pole p^ be 

G(s) ! G . ♦ ! G , ♦ G„ ♦ ... C3.U3) 

(s - Pi> (» - Pj) 

Then the angles of departure are given by 

“depart ' J l-‘a! * J • 0, 1. ... a - 1 

a = 1, 2, ... n (3.114) 

where X is a non- zero solution to 
a 


-n 


DET 


•(a>2) 


-(a*l) 


= 0 


G 


-(a*2) 


(=.a - 'a« 


‘^-(a.l) “-(a-2) 


-n 


(3.115) 


provided this is an equation in X^, and certain matrices obtained 
from the G^ by a procedure given in [17] have simple null structure. 


Theorem 3.10 is proved as follows. By making the substitution 
g » - in the characteristic equation, we have the following equation 
describing the root locus: 


DET [I + kG(s)] » 0. (3.116) 

This is equivalent to stating that there exists a non- zero m-vectorv(s) 
such that 


(I + kG(s))v(s) « 0. 


(3.117) 
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In the vicinity of a pole, where s - is near zero, we may expand v(s) in 
a power series 

v(s) a V^(s - p^) + v,(s - p^)" + ... (3.118) 

with V non- zero, 
o 

We now make use of the fact that if the all have simple null struc- 

1/ A 

ture, all loci departing from p^ depart in integer orders, i.e. as k 
where a is a positive integer. This is proved in Sastry and Desoer [17] 
for the case of asymptotes of root loci, and a similar argument applies for 
the angles of departure. So we may write 


Pi 


.1/a 


s s x^k , a = 1 . . . n 


(3.119) 

(3.120) 


for some constants x , 

a 

First, let a » n. Substituting (5.113), (3.118), (5.119), and (3.120) 
in (3.117), we get 

~n 


(V ♦ V 2 ♦ ... ) + i_ (i- G + G , 
0 1 n ^,n -n ^n-l -n+1 

x_ s s 

n 

* 0 . 


+ ... )(v^ v^2 + ... ) 


(3.121) 


Letting s-*-0 for the angles of departure, we require the constant term of 
(3.121) to be zero: 


+ -i- G 
n ■ 

^n 


n'^o 


(I 


— G )v 
n -n-' o 


0 . 


(3.122) 


Now let a = n - 1. This gives 



(V ♦ v.i ♦ ... ) ♦ G ♦ G * ••• )(V- ♦ V.s ♦ ... 

0 1 n-i ~n -n »n-i -n+i o i 


* 


and we require both the constant term to be zero 


(3.123) 


V * ■ G ,v ^ , G V, 
0 n-1 -n+1 o n-1 -n 1 


* n-1 ^-n+l^'^o n-1 ^-n'^1 

^n-1 \-l 


and the term of order to be zero 


(3.124) 


-i-r G V « 0. 
n-1 -n 0 


(3.12S) 


We may combine (3.124) and (3.125) into 


=-n n rv 

“J ['-1. 


(3.126) 


But must be non-zero, which can only happen if is such that 


* ’'nil') 


(3.127) 


There are two possibilities. Either there is no such in which case 

there can be no loci departing as or there are n-1 such 


specifically the n-1 solutions to 


n-1 . 

Vi ■ ■ ^ 


(3.128) 


where X solves 


DET 



0 


■ 0 . 


(3.129) 


".nj 

It is innnediatcly evident that the solutions to (3.128), when sub- 
stituted back in (3.120), will give rise to an (n - 1) -order Butter- 
worth pattern. Continuing the argument for a « n - 2, ... 2, 1, it is 
clear that the pattern of (3.122) and (3.127) will continue, giving rise 
to (3.115). This concludes the proof. 

For the angles of arrival, we have 


Theorem 3,11 Let the Laurent series expansion of G(s) at an n^^- 

order zero z. be 
1 

G(s) » (s - ♦ (s - * “* (3.150) 

and let the (again, see [17]) all have simple null structure. Then 
the angles of arrival are given by 


« arrival * ' I I'^al * - 3 * <>. 1. • • • 


a - 1 


a ^ 1, 2, ... n 


(3.131) 


where \ is a non-iero solution to 

3i 


DET 


a«2 


a-l 


a-2 


(G 


a 


^a'^ ^a-1 


V2 


0 

G. 


(3.132) 


provided this is an equation in X.* 

a 


The proof of Theorem 3.11 follows closely the proof of Theorem 3.10, 


and should be quite apparent 
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Solving the generalized eigenvalue problems (3.115) and (3.132) is by 
no means trivial, and from a computational standpoint our earlier results 
will often be preferable. The importance of Theorems 3.10 and 3.11 is 
the linkage they provide between all the results of this chapter. We now 
give an algorithm which shows dramatically the link between Theorems 3.10 
and 3.11 and Theorems 3.4 and 3.5. 


Lemma 3.1 Let the Laurent expansion of G(s) at a pole or zero q be 
1 ^ 1 


G(s) i— G^ ♦ - 

(s - q) ® Cs - q) 


^^1 ^ *** * ^c 


+ (s - * *** (3.133) 

and define the set of Toeplitz matrices 


i » 0, 1, ... (3.134) 


Now define the sequence 



« RANK - RANK i » 0, 1, ... (3.135) 

with RANK T ^ ■ 0. Then the structure indices of q as a pole are 
{k^, i « 0, 1, ... c • l) and the structure indices of q as a zero 
are {k^, i » c ♦ 1, c ♦ 2, ... 2c }. 


This result is due to Verghese [12] and Van Dooren ^ [13], and it 

can be used to compute the structure indices of G(s) without putting G(s) 
in Smith •MacMillan form. 

Now compare the matrices in (3.115) and (3.134). Except for the XI 
in (3.115), they are the same! Comparing Theorem 3.10 and Lemma 3.1 shows 


1 


that Theorems 3.4 and 3.10 are predicting the same orders for the Butter- 
worth patterns for the angles of departure. Similarly. Theorems 3.5 and 
3.11 predict the same orders for the angles of arrival. 

This agreement seems to indicate that the simple nuJl structure 
assumptions in Theorems 3.10 and 3.11 and the assumptions involving the 
principal minors of W(s) in Theorems 3.4 and 3.5 are equivalent. This 
would be a nice result, since it would confirm that the assumptions of 
Theorems 3.4 and 3.5 are indeed generic. Unfortunately, we have not been 
able to show this equivalence, although we suspect strongl<' that it exists. 

Theorems 3.10 and 3.11 can also be related to the iTurent series 
results of Secti<'n 3.4. If p^ is a first-order pole, only one branch may 
depart from it, and hence there can only be one non-zero solution to 
(3.115). For a first-order pole (5,115) becomes 

DET - Xjl] » 0 (3.136) 

and so Xj is the single non-zero eigenvalue of But TR G_^ is the sum 

of the eigenvalues of G^^. and if all of them except X^ are zero, then 
(3.114) yields 

® depart “ ^ ^’^1^ * ^-1^ 

which agrees with Theorem 3.6. 

If p^ is a simple n^^-order pole. (3.115) becomes 

DET [G „ - X„I] • 0 (3.138) 

-n n 

and again there can only be one non-zero solution. This will again be the 
single non-zero eigenvalue of G^^. and (3.114) now yields 
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« d.p.rt ■ K ““ 1-N.l 


♦ t j *0, I, ... n - I 


■ i ARC l-TR j ■ 0, 1, ... n - I (3.139) 


which agrees with Theorem 3.8. 

As this thesis was being written, we were made aware of concurrent 
and independent research into the natures of the angles of arrival and 
departure by Byrnes and Stevens [20]. In [20], Byrnes and Stevens derive 
the main result of Section 3.3.2, i.e. that loci depart from poles and 
arrive at zeros in Butterworth patterns whose orders are the MacMillan 
indices of G(s) at the pole or zero in question, lltey also show that the 
assumptions of Theorems 3.4 and 3.5 are in fact equivalent to simple null 
structure assumptions on the matrices obtained when G(s} is block- 
diagonalized, and hence are generic. However, they derive no explicit 
equations for the angles of arrival and departure, and they do not consider 
the approaches taken in Sections 3.4 and 3.5. 



CHAPTER IV 


Branch Points and Break Points 


4.1 Introduction 

In this chapter we consider the other two types of ^oints, besides 
poles and zeros, that are significant to the behavior of the multivariable 
root locus. These are branch points, which are associated with unusual 
behavior of the multivariable root locus and are not present in the SISO 
case, and break points, short for breakin points and breakout points, where 
a branch on the real axis suddenly breaks out into the complex portion of 
the s-plane, or the reverse. 

Branch points are perhaps the most startling new phenomenon encountered 
in generalizing the root locus from the SISO case to the multivariable case. 
They are associated with the unusual loops and swerves sometimes seen in 
multivariable root loci (see [5]). In particular, the branch points on the 
real axis are associated with the "turnaround" of root loci on the real 
axis depicted in Example 4.1 below. It will be shown in this chapter how 
branch points may be computed and how the "turnaround" behavior may be 
predicted. 

Break points are well known from the SISO root locus, but they have 
not been considered in the multivariable context. We show in this chapter 
how break points may be computed, and that the angles at which branches 
break into and out of the real axis are evenly spaced over 360®, as in the 


SISO case. 




— TrggJ IBKlitfii 'iLXgli ^'7' 
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4.2 Branch Points 

4.2.1 Computation of Branch Points 

Before the effects of branch points on the root loci can be ascer- 
tained, the branch points themselves must be computed. We now give two 
procedures for computing the branch points: one for the case of a two- 

input-two-output system (m » 2) , and one for systems with three or more 
inputs and outputs. 

Theorem 4.1 If m = 2, the branch points are given by the solutions to 
A(s) = (TR G(s))^ - 4 DET G(s) * 0 (4.1) 

and the gain g^ at a branch point s^ is given by 

go=iTRG(s^). (4.2) 

Theorem 4.1 follows immediately from the characteristic equation of a 
system with m * 2, which is 

g" - (TR G(s))g + DET G(s) = 0. (4.3) 

Solving this, we have 

g = l^TR G(s) ± V(TR G(S))‘ - 4 DET G(s)^ (4.4) 

and recalling (from Chapter II) that branch points are by definition points 
where the characteristic equation has a multiple root g^, the result follows. 

Theorem 4.2 The branch points of a root locus described by 4(g,s) = 0 
are given by the solutions to the simultaneous equations 

4(g,s) = 0 (4.5a) 

4(g,s) » 0. (4.5b) 


Remark The resultant (see Appendix) may be used to solve (4.5). 


We prove Theorem 4.2 by noting that if is a multiple root of 
» 0, we may write, for some s^, 

^(«.Sq) ■ (g - gQ)‘^(g.s^j) (4,6) 

and we have 

ig(g.Sg) » :j(g - ^ - So^^Tg^^.s^). (4.7) 

Setting g » gg makes both (4.6) and (4.7) zero. 

It should be noted that higher-order poles and zeros may be branch 
points (with * 0 or -«). We exempt this case from the results to 
follow. 


4.2.2 Effects of Branch Points on Root Loci 


The effect that a branch point can have on the form of the root locus 
is best illustrated by an example. 


Example 4.1 Plot the root locus for 

m m 

It is shovm in Chapter V (and also in [1], from which this example is 
taken) that the root locus is as drawn in Fig. 5.2. Note that the 
branch departing from the pole at -1 moves in the positive real 
direction, then abruptly turns 180* at the branch point at 1/24. 
Certainly this type of behavior is not characteristic of SISO root 
loci! 


As explained in Chapter II, this unusual behavior is observed at branch 
points because it is at these points that the root locus "jumps" from one 
sheet of its Riemann surface to another, which maintains continuity but 
allows for a sudden change in direction. This is discussed in more detail 
in [1] using the 180* phase contours of the algebraic function g(s); here 


we are more interested in describing this behavior than in accounting for 
it. 

It should also be noted that the root locus can "jump" from one 
Riemann surface to another at places where it crosses a branch cut . Branch 
cuts are the "seams" where different copies of the complex plane have been 
"stitched together" to form the Riemann surface. They are made between two 
branch points, or between a branch point and infinity, by a procedure 
described in [1]. The very unusual behavior of some of the root loci in 
[5] may be associated with branch cuts, but a detailed explanation of this 
behavior will require more research. 

The following argument, due to Postlethwaite [18], may be helpful in 
understanding why branch points produce the "turnaround" effect on root 
loci. Since '^'(g,s) = 0 for all (g,s) on the root locus, we have 


4 Hg,s) . i Hg.s) » ^ «(8.s)a| = 0 (4.8) 

and at a branch point (g^,s^) using Theorem 4.2, we have 

= 0 or (4.9a) 


ds I 

3 * 


0 . 


(4.9b) 


But (4.9b) implies that s^ is a stationary point of the root locus — a 
point where a branch turns around and doubles back on itself. 

We now show that it is possible to determine on which "side" of a branch 


point a branch of the root locus will approach, reach the branch point, 
turn around, and depart. This result is quite important since without it 
we can say little about the locations of loci on the real axis. The 


following theorem will be used in Chapter V: 


Theorem 4.5 Given a branch point on the real axis, the root locus 

will approach it, turn around, and depart from it on the left side 
(respectively on the right side) if 

SGNj-^-j— Is » s 1 * ^ (respectively -1). (4.10) 

L 3g 3s oj 

Remark Recall that we have g^ » iTR G(s^) at a branch point. 

We prove Theorem 4.3 as follows. In the vicinity of the branch point 

s„, define 
o 


6s * s - s^ 
o 


(4.11) 


and for a small perturbation jg in g write the Taylor series 

2 

6s = ~ 1 6g + J ^ I (6g)“ + ... (4.12) 

dg o dg o 

The first term is zero. Neglecting higher-order terras, one has 

SGN [6s] » SGN [~ ]. (4.13) 

dg*^ *0 

This illustrates the "turnaround" behavior -- regardless of the sign of 6g, 
which differentiates "approaching" and "departing," the sign of 6s does not 
change. Thus s is always on the same side of s^. 

From (4.8) we may write 


3 * 



3s 


and taking the derivative with respect to g, we get 


fz 


. / d li\3i/d 

d_s ^ L3s Vig 3g / ~ 5g \Ig as/J 

(Ilf 


dg 


and evaluating this expression at the branch point s^, we have 


1 - 
dg *0 



s * s^. 

c 


Using (4.13) and (4.16), we finally obtain 
SGN [6s] = -SGN 


9^4 ^ I 
2 s = s 

Lag 3s ®oJ 


(4.15) 


(4.16) 


(4.17) 


,3g 3s ■ 'O. 

and the theorem follows. 

There is an interesting corollary to this theorem in the special case 
m s 2: 


Corollary 4.1 If m « 2, (4.10) simplifies to 

SGN A(s^)j = 1 (respectively -1) (4.18) 

oS O 

if the loci are on the right side (respectively left side) of s^. 


This is not difficult to show. For m 
is (4.3), and (4.10) becomes 


2 the characteristic equation 


SGN 


3g^ as 



DET G(s)) 


(-TR G(s)g 



(4.19) 


But recalling (4.2) from Theorem 4.1, this becomes 
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SON 


l-3=N[-2((i- 
Lag- 3s ® oj L V 3s 


TR G(s))(1TR G(s)) 


- ^ OET G(s) 


- 4 DET G(s) 


)'* -“o] 


-SGN 


(TR G(s))^ 
(4.20) 


and the result follows. Thus the loci will be on the side of the branch 
point for which A(s) is positive -- a fact that we will interpret in 
Chapter V. 


4.3 Break Points 

4.3.1 Computation of Break Points 

We now give a procedure for computing the break points. Recall (again 
from Chapter II) that break points are by definition points where the 
characteristic equation has a multiple root s^. 

Theorem 4.4 The break points are given by the solutions to the simul- 
taneous equations 

*(g.s) = 0 (4.21) 

*(g,s) = 0 (4.22) 

Remark Again, we may use the resultant (see Appendix) to solve (4.21) 
and (4.22) by rewriting <t’(g,s) as 

ng,s) « y(s,g) ^ Bj^(g)s'^ + ®n- 1 ^ 8 )s"’^ + ... ♦ B^(g)-0. (4.23) 

Now the resultant yields an equation in g. We may then obtain the 
break point s^ by solving (4.22). 

Theorem 4.4 can be proved by repeating the proof of Theorem 4.2 with g 
and s interchanged, but the following proof is more insightful. Define 


(4.24) 


for (g,s) near (ggjS^). Since 'J'(g.s) * 0 for all (g,s) of the root locus, 
we have 

•SngQ.So) = H (80**0^'^® ^ li (So*V'^« * ° 

where we have taken 

6s s s - (4,26a) 

o 

<5g * g “ gjj (4.26b) 

sufficiently small to neglect higher-order terms. If we now consider a 
locus point on the real axis, all quantities in (4.25) will be real except 

6s = RE [6s] + j IM [6s]. (4.27) 

Substituting (4.27) in (4.25) and equating real and imaginary parts to zero, 
we get 

IM [«s] H Cg^.s,,) • 0 
«E [«sl II (g„,s„) . Sg II (g^,s^) . 0. 

If we now make s^ a break point, one has IM [6s] ^ 0 for the branch 
breaking in or out. This and (4.28) prove the theorem. 

4.3.2 Breakin and Breakout Angles 

We now show that the angles of loci breaking in or out are the same as 
in the SISO case. 


(4.28) 

(4.29) 


Theorem 4.5 If several branches of the root locus are approaching and 
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leaving a point on the real axis, their angles are evenly distributed 
over 360” , and branches approaching and leaving the break point are 
interleaved (i.e. they alternate). 

Remark Note that branches on the real axis must be included — not 
just the branches breaking in or out. 


To prove Theorem 4.S, we first show that if k branches are approaching 
a point s^ and k branches are leaving s^, then the first k - 1 derivatives 
of the algebraic function g(s) are zero at s^. We have 


d® , » ^ ^ dg I rt 


34 | 


3g 


(4.30) 


34 34 

If is a break point, (8 q»Sq) is zero and non-zero, so 
that 


d£ 

ds 


= 0 . 


(4.31) 


0 

If k different loci are all passing through s^, s^ is a root of multipli- 
city k of the characteristic equation. This means that we have 


(i) 

^-rry <f(go.s^j) = 0, i = 0, 1, ... k - 1. (4.32) 

i C * ^ 


By repeatedly taking derivatives with respect to s of (4.30) and using 
(4.31) and (4.32), it can be shown that 


- — p-r? «0, i =1,2, ... k - 1. 
ds^^^ 

Now, in the vicinity of the break point s^ define 


(4.33) 


6s = s - s^ 


(4.34) 


and for a small perturbation 6g in g write the Taylor series 
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*^8 • If Is ^ M U * ••• 

o ds o 


(4.3S) 


Of course the first k - 1 terms of (4.35) are zero, from (4.33). Taking 
the first non-zero term and neglecting the higher-order terms, we have 




(4.36) 


Also, by repeatedly taking derivatives of (4.30) we may obtain an expres- 

I 

Sion for — rt* I in terms of various partial derivatives of 4(g,s) 
ds'-'^^ o 

evaluated at (g ,5^), and these will all be real for s^ on the real axis, 
'•'o’ 0 — — 0 

d^^ 2 I 

Hence — | is a real number, and from (4.36) 
ds^*‘^ ®o 


k ARC [6s] • ARC [6g] - ARC L 


« 0‘ 

» 


or 180*. (4.37) 


We may now write 


ARC [6s] * 7 (9 ♦ n360“), n » 0, 1, ... k - 1 


(4.38) 


where 6 is 0® or 180®, depending on the sign of 6g. To find the angles for 
branches departing from s^, take 6g > 0, and to find the angles for branches 
arriving at s^, take 6g < 0. These form two k^*^-order Butterworth patterns, 
with principal angles of 0* and 180®. This proves the theorem. 
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CHAPTER V 

Root Loci on the Real Axis 



S.l Introduction 

It is generally very difficult to plot the root locus precisely for 
finite gains. Exact analytical expressions for .the various branches are 
usually difficult or impossible to obtain, and attempts to discern the 
locus by actually plotting the closed-loop poles for various values of k 
tend to be onerous at best. These difficulties hold even in the SISO case; 
they are considerably greater in the multivariable case. 

There is, however, one part of the root locus that can be plotted 
exactly with relative ease — the portion that lies on the real axis. The 
form of the locus on the real axis is of course known exactly, and, in 
addition, the number of branches of the root locus on the real axis can 
change only at a finite number of points. Thus a relatively small amount 
of work may yield an exact plot of a sizable portion of the root locus; 
indeed in some cases all of it (see Example 5.2). Knowledge of asymptotes 
and angles of arrival and departure is often sufficient to sketch the rest 
of it. 

In the SISO case the rule for the location of root loci on the real 
axis is very simple (see Chapter II). This is because only one branch of 
the root locus can lie on the real axis at any given point. However, in 
the multivariable case as many as m branches can lie on the real axis at 
a given point. Thus the problem is not one of determining whether a branch 
is present, but one of determining how many branches are present. The 
“turn-around" behavior associated with branch points (not present in the 



SISO case) makes things even more complex. 

The first observation is that, unlike in the SISO case, knowledge of 
pole and zero locations alone is not sufficient for determining the nundier 
of loci on the real axis. The following example makes this clear: 

Example 5.1 Plot the root loci for 


0 

s+z 


0 1 


n » - 2 

0 ■ ■■ r 

and G^Cs) « 

o 

II 

S • 1 

• « 


L * - 'J 


Since each of these represents two decoupled SISO systems, we may 
immediately plot the root loci, which ai*e given in Fig. 5.1. Note 
that although G^(s) and G2(s) have their poles and zeros at the same 

locations, the number of loci on the real axis between -1 and 1 are 
different . 

Despite this difficulty, some equations for the number of branches of 
the root locus on the real axis at a given point may be found. Also, these 
equations are not too complicated to be useful. We consider first the case 
m « 2, then the general case, and finally the case when G(s) is symmetric. 

S.2 The Case of Two- Input-Two-Output Systems 

The following theorem provides a step*by>step procedure for deter- 
mining the number of root loci on the real axis by solving three simple 
inequalities for bounds on s. 

Theorem 5.1 For a system with m ■ 2, let s be any point on the real 
axis. Define 

i(s) ^ (TR G(s))^ - 4 DET G(s). (S.l) 

Then we have: 

(i) If DET G(s) < 0 . exactly one branch lies on the real axis at s; 





(ii) If PET C(s) > 0 . two or zero branches lie on the real ffxis at s: 

(a) If ^(t) < 0 , zero branches lie on the real axis at s; 

(b) If i(s) > 0 and TR GCs) > 0 . zero branches lie on the real 
axis at s; 

(c) If i(s) > 0 and TR G(s) < 0 . exactly two branches lie on 
the real axis at s. 

To prove Theorem 5.1, we will follow the order in which it is stated. 
This will minimize confusion and also minimize the amount of work needed, 
since DET G(s) > 0 is necessary to have <i(s) < 0, and since TR G(s) is 
sometimes not needed. 

Ne start by observing that the characteristic equation with m « 2 is 

u(g,s) • g“ - (TR G($))g ♦ DET G(s) « 0. (5.2) 

Now let s be real and vary over the entire real axis. For a given s^, the 
number of branches on the real axis at s^ is equal to the number of negative 
real roots of (5.2) with s « s^ (recall that a negative real g corresponds 
to a positive real gain k). Since the roots of (5.2) are given by 

g • i (tR G(s) 2 VSTsl) (5.3) 

we merely investigate how many negative real values of g we get for various 
values of <l(s), TR G(s), and DET G(s). Thus if d(s) is negative the two 
values of g will be complex, and there will be zero branches on the real 
axis at s. But d(s) can only be negative if DET Q(s) is positive. If 
DE*r G(s) is negative, then 

iTR G(s)i < V(TR G(s))^ - 4 DET G(s) (5.4) 

and the two values of g are real and of opposite sign. Hence there is 


exactly one branch on the real axis at s. The other rules follow similarly. 

We now make some remarks that will hopefully aid in understanding some 
of the features of this theorem: 


1. Note that the number of loci on the real axis changes by one 
whenever DET G(s) changes sign. This makes sense, since loci 
start at poles and end at zeros and since DET G(s) changes sign 
at a (first-order) pole or zero. 

2. The number of loci on the real axis may change by two whenever 
A(s) changes sign. But the points at which A(s) will change 
sign are the branch points (recall from Theorem 4.1 that A(s) is 
zero at a branch point) , and the "tum-around" behavior of a 
root locus branch at a branch point would indeed make the 
number of loci change by two. In fact. Corollary 4.1 correctly 
predicted that the "side" of a branch point from which a branch 
would approach, turn around, and depart is the side for which 
A(s) is positive! 


3. The number of loci on the real axis will occasionally change by 
two at points where TR G(s) changes sign. This happens when 
there is a double pole or zero with both branches departing or 
arriving on the same side. For example, consider 


G(s) * 


s - 1 
s - 1 


s + 1 J 


Clearly there will be two branches both departing the pole at -1 
in the positive direction. We have 






■ ■«--- ■ 
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TR G(s) 


2s - 3 


s ♦ 1 

which does indeed change sign at -1, 


Note that if 


G(s) = 


^(s) 




we have 


TR G(s) = ^ + n 22 (s)) 

’ d(s)^ (nijCs)n22(s) - nj2(s)n2^ (s)) 

' d(s)2 (C^llCs) - n22(s))^ + '»ni2(s)n2i(s)) 


(S.S) 


CS.6a) 

(S.6b) 

(5.6c) 


and since d(s)2 is always positive (except at poles) we may neglect it in 
solving the inequalities for bounds on s. 

«e end this section with an example (taken from [!)) to illustrate how 
Theorem S.l may be implemented. 


cAainpie Plot 


*»*4w AWi. XUWUS 




G(s) * — ^ 

(s + l)(s + 2) 


s - 1 


-6 


s - 2 


i?Z““!:a^^‘i^’fi"XdTpo^^rar!r»7s 

zeros. Using (5. 6b) one has ^ t 1 and -2, and no finite 


2 

DET G(s) a g ■*• 3s » 2 

((s + l)(s + 2))^ 

so that DET G(s) <0 for -2 < s < 
branch on the real axis for -2 < s 


-1 and consequently there is one 
^ • 1 • 
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Using (5.6c) we have 

1 - 24s 

A(s) T 

((s + DCs ♦ 2))^ 

so that A(s) <0 for s > 1/24 and therefore there are no branches on 
the real axis for s > 1/24. 

Finally, using (5. 6a) we have 
' “ (s ♦'l)(s^+ 2) 

and restricting ourselves to the range of s not already considered we 
have TR G(s) <0 for s > -1 or s < -2, and this implies that there are 
two branches on the real axis everywhere else. 

The root locus branches on the real axis are plotted in Fig. 5.2. But 
since we have two poles, no finite zeros, and two asymptotes, this is 
the complete root locus! 


5.3 The General Case 

The general case when m > 2 is much more complicated than the case 
when m * 2. However, after evaluating a few quantities, we may use the 
following theorem to determine immediately the exact number of loci every- 
where on the real axis. 


Theorem 5.2 Assume that all higher-order poles and zeros on the real 
axis are "simple" ("simple" poles and zeros were defined in Chapter 
III) . Then the exact number N of loci on the real axis at a point s^ 
is given by 

/ number of \ 


N 


zeros Zj^ of 
odd order to 
right of So 


f k. 

SGN (s - z^) ^ TR 


1 / number of \ 

G"^(s)l^ . 2 j + (asymptotes j 


poles Pi of 
odd order to 
right of So 


+ 2 


branch points 
bi to right 
of s„ 


k. 

SGN (s - pp ^ TR G(s)l^ 




.,1 


(5.7) 



where the are the orders of the poles and zeros, and where the 

summations are taken over zeros and poles of odd order and branch 

points on the real axis and to the right of s^. 

It should be noted that (S.7) is not nearly as complex as it might 
appear at first glance. In order to apply Theorem 5.2, we need only 
evaluate the sign of a quantity at each pole, zero, or branch point. 

Once this has been done, we simply add up the different contributing terms 
of (5.7) for each stretch of the real axis between any two of the three 
types of points to determine the number of loci on that stretch of the 
real axis. 

•ve prove Theorem 5.2 by a "conservation of loci" argument: each locus 

must start somewhere, end somewhere, and be continuous in between. We start 
with just poles and zeros, and then consider break points, asymptotes on 
the real axis and branch points, adding them in as we proceed. 

We claim first that if there are only first-order poles and zeros on 
the real axis, the number of loci N on the real axis at a point s^ is 
given by 

j, ^ .niimber of poles to right of Sq •> 

* '■with a branch departing at 180"'^ 

.number of zeros to right of Sqn 
^ with a branch arriving at 180®^ 

number of zeros to right of Sq^ 

* 'with a branch arriving at 0" ^ 

.number of poles to right of Sq. re «■» 

^with a branch departing at 0® •* * J 

This is easy to see, since the first two terms give the number of 
branches moving in the negative real direction and the last two terms give 
the number of loci moving in the positive real direction at s^. (Recall that 


the angle of arrival is the direction in which the locus is moving when 
it reaches a zero.) 

We now extend this to higher-order poles and zeros that are simple. 
Recall that a simple pole or zero will have loci departing from it or 
arriving at it in a single Butterworth pattern. By symmetry, we see 
immediately that a simple pole or zero of even order can have no effect 
on the number of loci on the real axis, since such a pole or zero will 
either have loci departing or arriving at 0* and 180*, or no loci departing 
at either 0* or 180*. Either way, there can be no contribution to (5.8). 

On the other hand, a simple pole or zero of odd order must have exactly 
one locus departing or arriving at either 0* or 180*. The simplest way 
to determine the angle is to use Thef^rems 3.8 and 3.9, and upon substitution 
in (5.8) these yield the first and third terms of (5.7). Note that since we 
are only considering points on the real axis, all quantities will be real, 
and in (5.7) we may use the SGN function instead of the ARC function. 

Next, we introduce break points and asymptotes on the real axis. It 
is easy to see that break points have no effect on the number of loci on 
the real axis, since breakin points act like double poles with loci 
departing at 0* and 180* and breakout points act like double zeros with loci 
arriving at 0* and 180*, and neither of these makes any contribution to 
(5.8). However, asymptotes on the real axis are zeros at infinity with 
branches arriving at 0*, so according to (5.8) we must a^ them in. This 
yields the second term in (5.7). 

Finally, we must introduce branch points. Here we must use Theorem 
4.3, which stated that the "side'* from which a branch will approach a 
branch point b^, turn around, and depart from it is given by 



depending on whether the locus is on the left side (respectively the right 
side) of b^. If the locus is on the left side of b^, then b^^ is acting like 
a combination of a zero with a branch arriving at 0** and a pole with a branch 
departing at ISO”. In that case, according to (5.8), we should add two to 
the number of loci on the real axis if b^ is to the right. On the other 
hand, if the loci are on the positive side of b^^ everything is reversed 
and, according to (5.8), we should subtract two from the number of loci. 

We use Theorem 4.3, and this yields the final term in (5.7). We have now 
considered all possibilities, and Theorem 5.2 is proved. 

The following corollary is interesting, primarily because it is the 
closest we can come to generalizing the SI SO rule for loci on the real axis 
to the multivariable case. It may also be used as a check when applying 
Theorem 5.2, and may even provide sufficient information by itself for 
some applications. 

Corollary 5.1 Assume that all higher-order poles and zeros on the real 
axis are simple, and that there are no asymptotes on the real axis at 
+®. Then, counting multiplicities, at least one branch (in fact, an 
odd number of branches) of the root locus lies on the real axis at 
a given point s^ if there is an odd number of poles and zeros to the 

right of s^. 

Remarks (1) There is an odd number of poles and zeros to the right 

of 5^ if and only if DET G(s^) < 0 (this is proved below); 

(2) If there is an even number of poles and zeros to the 

right of s^, then there is an even number of branches on 

the real axis at s^. Unfortunately, zero is an even 

number. 


Corollary 5.1 follows almost immediately from (5.8). Making the 


obvious substitutions, write (5.8) as (x^ - X 2 ♦ Xj - x^) . The total 
nun^er of poles and zeros to the right of s^, counting multiplicities, is 
(x^ ♦ X 2 ♦ Xj ♦ x^) , and it is clear that (5.8) will be odd if and only if 
this quantity is odd, guaranteeing at least one branch on the real axis at 
s^. Recalling that break points have no effect on the number of loci on 
the real axis, and that branch points can only change the number of loci 
by an even number con^letes the proof. 

The first remark follows from the fact that DET G(s) changes sign at 
the poles and zeros of G(s) (counting multiplicities), and from the following 
argument which proves that DET G(s) is positive at infinity if there are no 
asymptotes there. Suppose that DET G(s) is negative at infinity. Then the 
product of the eigenvalues of G(s) at infinity is ne“?ative, and G(s) must 
have at least one negative real eigenvalue at infinity. But the charac- 
teristic equation for the root locus is 

A(g,s) » DET (gl - G(s)) » 0 (5.10) 

so if G(s) has a negative real eigenvalue at infinity there must be a 
branch of the root locus at infinity, contradicting the assumption that 
there are no asymptotes on the real axis at **». 

There is another way of showing that if DET G(s^) is negative there 

is at least one branch of the root locus on the real axis at s^. Consider 

o 

A(gjS^) to be a polynomial in g of degree m with constant term 
(-1)™ DET G(Sq). If m is odd and DET C(s^) is negative, then the constant 
term is positive. But for large, negative values of g, A(g,s^) will be 
negative. Since A(g,s^) is continuous, it must cross the negative g-axis 
somewhere, and hence A(g,s^) must have at least one negative root, which 
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implies that a branch of the root locus lies on the real axis at s^. If 
m is even the same argument applies, with the signs of (>1)"^ DET G(s^) 
and A(-»,s^) reversed. 

Considering d(g,s^) in this fashion also gives us an i;^per bound on 

the number of loci that may lie on the real axis at s^. Since ACg.s^) has 

degree m, it can have at most m negative real roots. Hence there can be 

at most m loci on the real axis at s . 

o 

S.4 The Case of Symmetric G(s) 

In this section we specialize to the case when G(s) is symmetric. 

Since G(s) is symmetric for reciprocal networks, this case does have some 
practical applications. Our final result depends on matrices obtained from 
G(s) by several transformations, so for clarity we will proceed with its 
derivation and then, having derived it, state the result as a theorem at 
the end of this section. This result is due to Levy [27]. 

The characteristic equation of the root locus is 

d(g,s) » DET [gl - G(s)] » 0 (5.11) 

so the number of branches of the root locus at a point s on the real axis is 
the number of negative real eigenvalues of G(s). However, if G(s) is 
symmetric then all of its eigenvalues are real , and we need only to 
determine how many of them are positive and how many of them are negative. 

In order to keep track of the number of negative eigenvalues of G(s), 
we will use the signature of G(s), which we now define. 


Definition Let M be a n'n-singular real symmetric matrix, and define 
m^ s the number of positive eigenvalues of M (S.12a) 
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the number of negative eigenvalues of M. 


(5.12b) 


Then the signature a (M) of M is defined as 


a(M) « - m_. 

Remark Since M is non-singular, we have 


(5.13) 


m^ + m ■ m 


(5.14) 


where m is the size of M. Therefore we may determine m and m from 
o(M). 


The signature of a matrix is a useful concept in the present context 
because it is invariant under congruency transformations. Thus if L is a 
non-singular real matrix and we have 


P » LML 


(5.15) 


then o(P) « a(M). We will use this property several times in this section. 
Now write the left matrix fraction description of G(s) 


G(s) . D’^(s)N(s) 


(5.16) 


where D(s) and N(s) are left coprime polynomial matrices. The poles of G(s) 
are the zeros of DET D(s) , and the zeros of G(s) are the zeros of DET N(s). 
Since the product of the eigenvalues g^^Cs) of G(s) is given by 


• “ST • iHB} 


(5.17) 


the eigenvalues gj^(s) can only change sign at the real poles and zeros of 
G(s). 

For all points on the real axis that are not poles and zeros of G(s), 
we have 


D(s)G(s)d”‘‘(s) . N(s)D^(s) - P(s) 


(5.18) 



is a congruency transformation of G(s), and hence 

i. ■- 

‘JCG(S)) - o(N(s)d'^(s)) • a(P(s)) (5.19) 

Since: (1) the number of loci on the real axis at s is the number of 

4 negative real eigenvalues of G(s); (2) the number of negative rual eigen- 

; values of G(s) may be determined from o(G(s)) « <7(P(s)); and (3) cr(P(s}) • 

o(G(s)) can only change at a pole or zero of G(s), we now investigate 

how o(P(s}) changes near a pole or zero s^. A procedure for doing this 
follows. 

We may write 

P(s) * Pq + P^Cs - Sq) + ... + Pj(s - Sq)*^ (5.20) 

where the P^^ are real and symmetric and where P^ is singular. Hence there 
exists a real non-singular matrix such that 


o o 


L 






0 Oj 


(5.21) 


where G^ is real, symmetric, and non-singular. We may then define 



Fg o1 

d 

■A, B,' 

Q(s) « T P(s)tJ • 

0 

* E 


.0 0, 

i«l 

cJ 


(s - s^)^ 


(5.22) 


and Q(s) is congruent to P(s) . 

We now zero out by using another congruency transformation. Define 


''l - =o‘»l 


5 » s - s. 


(5.23) 

(5.24) 


I 0 I 

T Q(S) 

-vjs I 0 


(5.25) 


and note that we may write 


0 A. 0 

o Is 

S ♦ ... 

0 0 0 c. 


(5.26) 


where is real and symmetric. If also has full rank we may halt this 
procedure, since in this case (5.26) has the desired form given in (5.31) 
below. 

If does not have full rank, this procedure must be repeated 
starting with (5.21), and with taking the place of P^. That is, we 
write 




Gj 0 
0 0 


(5.27) 


where is a real non-singular matrix and where is real, symmetric, and 
non-singular. Now define the congruency transformation 


I 0 


I 0 


(5.28) 


and write 


G^ 0 A I 0 
® ♦-i-i 

0 o\ 

I 0 0 


d FaP) b«> 
cf> 


Next, zero out using a congruency transformation of the form 


(5.29) 
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T(s). ^ ° S(s) 

f-VjS^ - vjs) I 


I - Vjs) 


(S. 30 ) 


and if the resulting does not have full rank, run through the entire 

m 

procedure yet another tine. The procedure will terminate when we have 
obtained a polynonial matrix of the form 


G. 0 I A ' 0 


0 0 




I Gj 0 
I 0 0 


I G2 0 
I 0 0 


I 0 oJ L 


(5.31) 


(i+1) 

where » DIAG (A^ •* ,G^J , and the G^^ all have full rank. 

Nowwemay investigate how a(P(s)) changes near s^. Since Z(s) was 
obtained from P(s) by a series of congruency transformations, we have 


<^(^(s)) * <’(P(s)). In the vicinity of s^ Z(s) nay be approximated by 


5 U*1 

Z(s) a OIAG (G^,, Gj5, G2i , \ 


(5.32) 


and the eigenvalues of Z(s) are the eigenvalues of G^, G^s, ... and 

Now consider what happens to the signs of the eigenvalues of G^s^ if 
I changes sign from positive to negative. If i is even there will be no 
changes in the signs, but if i is odd all of the positive eigenvalues will 
become negative, and vice*versa. Then o(G^S^) will change sign, and so 
the change in c(G^5*') will be -2o(Gj^). (Note that if s is positive, 
o(G^ih ■ <^(Gj^).) It follows immediately that the change in ^(Z(s)) will 


9S 


be 

Ao(2(s)) • -2 iE o(GJ. (S.33) 

i odd ^ 

Prom (5.13) and (5.14) we may write 

mJZ(s)) • l(m - a(2(s))) (5.34) 

where m^(Z(s)) lathe number of negative eigenvalues of I(s) and m is the 
size of Z(s). This in^lies that the change in the number of negative real 
eigenvalues will be 

Am (Z(s)) « -iAo(Z(s)) « 2 o(G.). (5.35) 

i odd ^ 

Now let s vary along the real axis from to - ” , and assume G(») is 
positive definite (this is equivalent to assuming that there are no 
asymptotes on the real axis at *»). For each pole or zero s^^ on the real 
axis, we may compute a set of matrices G^^^ by using the above procedure. 
Then, recalling that a(Z(s)} • i7(P(s)) ■ a(G(s)) and that the number of 
branches on the real axis at s is the number of negative real eigenvalues 
of G(s) , we have proven 


Theorem 5.3 Assume that there are no asymptotes on the real axis at 
and that G(s) is symmetric. For each pole or zero s^ on the real 

axis confute the matrices cj^^, using the procedure described above. 

Then the number of branches N on the real axis at s is given by 


N. Z 

all poles and 
zeros s^ to 
right of s 


^ <j(g5^^). 

j odd ^ 


(5.36) 


There is an interesting observation that may be made on the procedure 


for generating the G^. Consider the set of Toeplitz matrices 


T(i) 



, i > Oj 1, ... d. 


(5.37) 


Recall from Lemma 3.1 that the structure indices of P(s) at s^ are given by 


» RANK T(i) - RANK T(i - 1) , i = 0, 1, . . . d. 


(5.38) 

However, it may also be shown that the congruency transformations used 
to generate the may be applied to the T(i) , yielding matrices of the form 


Go 0 ! I 

o * n 

I I ••• 

0 0 • i 

A, T 0 I G T 0 I 

-i-l 4-»+ 1 

I G, 0| I 1 ... 

0 I ^ f I I 

i ° i*ij ° i 

I G- 0 1 0 1 . . . 

0 I ; 0 ' 

0 0 I 0 0 




so that we have 


RANK G^ « RANK T(i) - RANK T(i - 1) , i « 0, 1, . . . d. (5.39) 

Therefore the pole* zero structure at s^ is determined by the ranks of 
the Gj^, while the number of branches of the root locus at s^ is determined 
by the signatures of the G^. 
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CHAPTER VI 


Asymptotic Behavior of Root Loci 
6.1 Introduction 

In this chapter we discuss the behavior of the branches of multi- 
variable root loci that tend to infinity as the feedback gain k approaches 
infinity. We have already examined the behavior of the branches that 
approach the finite zeros of G(s) as (i.e. the angles of arrival), and 
we now examine the behavior of the branches that do not approach finite 
zeros. These branches approach asymptotes as k^**, and it should be evident 
that knowledge of these asymptotes would be a considerable aid in plotting 
the root locus. 

An asymptote is characterized by its angle and by its order, which 
gives the "velocity" at which the locus tends to infinity. Asymptotes 
start at points called pivots, and unlike the SISO case, pivots for multi- 
variable root loci may be complex (see Exasq>le 6.1). 

We consider first the case where all asymptotes are of first order, 
since this case holds in general. The Newton polygon technique is used 
to show how the first-order asymptotes may be obtained from the eigenvalues 
of the first Markov parameter of G(s). We then give equations for the 
first-order asynptotes and pivots based on the characteristic equation, 
and give an exaaple with first-order asymptotes, complex pivots, and 
branch points, breakpoints, and loci on the real axis. 

The Newton polygon technique on which this thesis is based seems to 
be sonewhat inapplicable to the case of higher-order asymptotes. Also, a 
considerable amount of work has already been done on this subject, using 
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l 




I 


other approaches. For the sake of completeness, we conclude this chapter 
with a summary of the results of Shaked and Kouvaritakis [15] and Sastry 
and Desoer [17] on this subject, and we also note some other results. 


6.2 First-Order Asymptotes and Pivots 

In this section we examine the case in which all asymptotes are of 
first order. Since this is true in general, the results of this section 
are usually sufficient for computing all asymptotes and pivots. Ne show 
first how the asymptotes may be obtained from the first Markov parameter, 
and then how both the asymptotes and their pivots may be obtained fr<m the 
characteristic equation. 


Theorem 6.1 Let the Laurent series expansion of G(s) at infinity be 

G(s) « -i-Gj ♦ ~G 2 ♦ ... (6.1) 

s s 

Then, if and only if G^ is non-singular, all asymptotes are of first 
order with angles 


^asymptote “ i “ 1 ••• ® 

where the are the eigenvalues of G^. 


( 6 . 2 ) 


Remark If the system is described by state-space matrices (A,B,C), we 
have G, « CB. 


Theorem 6.1 is not a new result, but the following proof which is 
based on the Newton polygon technique gives an interesting picture of what 
is happening. Since we are interested in the behavior of s as S'***, sub- 
stitute s • 1/z in (6.1). This gives 

* G(s) « zGj ♦ z^Gj * ... * z(6j ♦ zGj +...)■ zG(z) (6.3) 

where 6(z) is a matrix power series in z. Substituting (6.3) in the 


] 


I 

i 

■ ■) 

i 

I 

j 

I 


I 

( 


i 




characteristic equation, we have 


A(g,s) - DET (gl - G(s)l 


g" 


z TR G(t)g®“^ 


♦ z 


E /principal minors \.m-2 
\of order 2 of 6(z)/* 


- ♦ 2®(-l)“ OET 6(2) « 0. 


(6.4) 


Since G(z) is a matrix power series, all of the coefficients in (6.4) are 
power series in z. The Newton polygon for (6.4) is drawn in Fig. 6.1, and 
it is clear that all of the points will lie on or above the line extending 
from the point (0,m) to the point (m,0). Since the coefficient of g® is one, 
the point (m,0) is definitely part of the polygon, and if and only if the 
constant term of DET G(z) is non- zero, the point (0,m) will also be part 
of the polygon. In this case, the Newton polygon will be the single line 
drawn in Fig. 6.1. 

Now note that the constant terms of DET G(z), TR G(z), etc. are just 
DET Gj, TR Gj^, etc., since 2(0) « Gj^ and since the constant terms are 
obtained by setting z equal to zero. Thus the Newton polygon will be as 
in Fig. 6.1 if and only if DET Gj^ is non-zero, i.e. iff Gi is non-singular . 

We then have 


g s cz for (g, 2 )-*-( 0 , 0 ) 
or equivalently, 

s a -ck for (s,k)-*‘(*,*) 


(6.S) 


( 6 . 6 ) 




where c solves 



c- - TR . (-»" BBT 0, 

« DET [cl - Gj] • 0. (6.7) 

This shows that c is an eigenvalue of G^. Applying the definition of angle 
to (6.6) concludes the proof. 

We now show how the first-order asymptotes and pivots may be obtained 
directly from the coefficients A^^(s) of the characteristic equation 
♦(g*s) ■ 0. 


Theorem 6.2 Given the equation 


m-1 


Hg,s) * Ajjj(s)g*“ ♦ A^_j(s)g"‘ * + ...+ A^(s) - 0 

we can write A^^(s) as 

a. (a.-l) 

A.(s) * g.s ♦ Y-s + (lower -order terms). 


1 1 
X ^0, 1, ... m 


and define 

M^{i:aj^-i>.aj-j, j«0, 1, ...m> 
N ^ {i; i + 1 C M}. 

Then the first-order asymptotes are given by 
Sj^ s p^ ♦ Cj^k, i > 1 ... m 


where the c^^ solve 


L(5.(-D^c{) 

jCM ^ ^ 


( 6 . 8 ) 


(6.9) 

(6.10a) 

(6.10b) 

( 6 . 11 ) 

( 6 . 12 ) 


and the pivots p^^ are given by 


Pi 


(6.13) 




Remark The angles of the asymptotes are given by ARG . 

The proof of Theorem 6.2 is straightforward. The first step Is to 
assume that 


s a p ♦ ck, (s,k)-^(«,“) (6.14) 

since this is the behavior of interest. Solving for k, we get 

k a . (6.15) 

Now substitute g « -1/k in (6.8) and multiply by k"*. This gives 

Ajj(s)k® - Aj(s)k®‘^ ♦ ♦ (-l)®Ajj^(s) » 0. (6.16) 

Substituting (6.15) in (6.16) and multiplying by c® gives 

Aq(s) (s - p)® - Aj (s) (s - p)®“^c ♦ . . . + a 0. (6.17) 


The degree of the i^** term on the left side of (6.17) is (a^^ + m - i), where 
the exponents a^ are defined in (6.9). The terms of highest degree are 
the terms, where j C M and where M is defined in (6.10a). Since the 
left side of (6.17) is asymptotically equal to zero, the sum of the 
coefficients of these terms must be zero. This yields (6.12). The sum 
of the coefficients of the terms of degree one less than the highest degree 
must also be zero. This gives 

V J (- 1) . 8 J (-1) J (. - J ) (-P)ci ) . E J 8 j (-1) JcJ ) . 0 


(6.18) 
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and solving for p yields (6.13). This concludes the proof. 


6.3 An Illustrative Example 

In this section we give an example that illustrates the computation 
of several features of the multivariable root locus. 


Example 6.1 Draw the root locus of 

» « 

, s 2s '«■ 6 

G(s) i—j 

(s ♦ 2)^ -(s ♦ 2) s ♦ 2 

Poles and Zeros 

We have 

2 2 

♦Cg>s) • (s + 4s + 4)g - (2s ♦ 2)g ♦ 3 » 0 

so there is a double pole at >2 and no finite zeros. Hence there are 
two asymptotes. 


(i) By inspection, we have 


1 2 


‘ -J 

which has full rank and eigenvalues 1 ± /? j. Hence, from Theorem 
6.1 there are two first-order asyn^totes with angles 

ARC [-1±/I j] » 12S.3*, 234.7*. 

(ii) By examining the coefficients of ^(g,s), we have 

M » (0, 1, 2), N - 4>. 

Using Theorem 6.2, (6.12) becomes 
2 

c ♦ 2c ♦ 3 • 0 

which has the solutions -1 ± j. Therefore the angles of the 
asymptotes are 

ARC (-1 ± /I j] • 12S.3*, 234.7*. 


2(2 - 1)(-1 ♦ j) ♦ 3(2 - 0) 


lOS 


which have the solution (Sg*s^) * (-1,-3). 

Since is real and negative, -3 is a break point. We know from 

Theorem 4.5 that the angle of the branch breaking in or out will be 
90*. 

We now have enough information to plot the entire root locus. The 
root locus plot is given in Fig. 6.2. 

6.4 Higher-Order Asymptotes 

We now review briefly the results of Sastry and Desoer [17] and Shaded 
and Kouvaritakis [15] on the subject of higher-order asymptotes. Proofs are 
omitted, since they do not employ the main methodology of this thesis. 

In Section 3.5 we adapted the procedure used by Sastry and Desoer in [17], 
and obtained equations for the angles of arrival and departure which involved 
solving a generalized eigenvalue problem in a Toeplitz matrix. We now state 
Sastry and Desoer *s original result, which dealt with asymptotes of root 
loci: 


G(s) ■ — G2 ♦ ... 


Theorem 6.3 Let the Laurent series expansion of G(s) at infinity be 

(6.19) 

i'*'. 

n * n 

n n * n n 

j ■ 0, 1, ... n - 1; n * 1, 2, ... 


s s 

(i) Then the angles of the n'’“-order asymptotes (if any) are 
^asymptote. 


( 6 . 20 ) 


where solves the generalized eigenvalue problem 


DET 






(6.2i) 


! 



provided this is a polynomial equation in 


(ii) If certain matrices obtained from the and defined below in 

(6.26) have simple null structure, then all of the asymptotes 
are integer>ordered, and hence are given*^ (6.20) and (6.21). 

If this assunq[>tion does not hold, then it is possible to have 
fractional-ordered asymptotes. 


(iii) 


The pivots c for the n^^-order asysqptotes may be found by 
solving " 


r 


0 


on 


DET 






'n-1 


• c G ) (G - X I) 
n n-' 


n 


0 

G 


IJ 


( 6 . 22 ) 


These results are of course proved in [17], although since virtually 
the same procedure was used to prove Theorem 3.10, it should not be diffi* 
cult to see where these results come from. 

In [17], Sastry and Desoer interpret these generalized eigenvalue 
problems as finding the eigenvalues of restricted linear maps. They also 
consider some ways of sioqilifying these problems, based on the Toeplitz 
structure, to facilitate their solution. The interested reader is referred 
to [17]. 

Comparing Theorems 6.3 and 3.10, and recalling how Lemma 3.1 linked 
Theorem 3.10 to the Smith-MacMillan form of G(s), the question arises as 
to whether there is any relation between the orders of asysqitotes and the 
Smith'NacNillan form. Verghese and Kailath [22] have pointed out that 
there is indeed a relation: the order; of the asymptotes are the orders of 
the Smith-MacMillan zeros of G(s) at infinity. The importance of this 
result is that the asysqitotes may be regarded as branches arriving at 
infinite zeros, so that the angles of the asynptotes are merely the ''augles 
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of arrival" at infinir« zeros. In earlier literature on the root locus 
(e.g. [15]), infinite zeros were considered to be fictitious objects. But 
they are in fact perfectly well defined from the Smith-MacMillan form at 
infinity. (For a discussion of the Smith-MacKilian form at infinity, see 
[23], p. 449.) 

The approsch taken by Shaked and Kouvaritakis [15] is completely 
different, and while their results are more thorough, they are also more 
arduous c«q)utationally. In [15] Shaked and Kouvaritakis approach the 
(troblem from a state-space perspective, and interpret their results in terms 
of mappir.^5 between spaces defined by the ranks and nullspaces of the Markov 
parameters We summarize their main results in the following theorem: 


Theorem 6.4 Let the Laurent series expansion of C(s) at infinity be 


G(s) ■ — G. ♦ G 2 ♦ ... 
s * s^ ^ 


and define the projected Markov parameters G . , i 
following sequrace: ^ 


G, 


hJ Taj ol 

[0 0] [ n ^ 

G^^^ * ^i^i 1 * * * ^l^i^l^l^2 * * * 


(6.23) 
. V using the 

(6.24 J 

(6.25) 

(6.26) 


where (6.25) is the spectral decosqposition of G^, exhibiting its Jordan 
form: (6.26) is the projection of «i.i onto the nullspaces of ... 

and we have assumed that all of the have simple null structure . The 
sequence terminates at i ■ v when G^ has full rank. 


The angles of the i^^-order asymptotes are then given by 

n360* 


i * T ‘‘i.j' * 


,n *0, 1, ... i— 1: 


j • 1. 2. 


*‘^i" ^ 2. ... V (6.27) 
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where X is the j non- zero eigenvalue of G. and d. is the rank 
deficiency of (d^ ■ m). 

Shaked and Kouvaritakis also give an alternative approach that can be 
used when the simple null structure assumption on the G^ does not hold. 
Kouvaritakis [24] discusses the use of a constant gain pre-compensator to 
make the G^ have simple null st.i ' ore, which reduces the orders of the 
asymptotes and thus improves gain margins. In [25] Kouvaritakis extends 
these results to the case of non-proper systems, and in [26] he applies 
these results to the optimal root loci associated with the linear quadratic 
regulator problem. 
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CHAPTER VII 
Other Results 

7.1 Introduction 

In this chapter we discuss briefly several other results on multi- 
variable root loci that have recently appeared in the literature. Although 
these results are rather minor, we feel that they are sufficiently 
interesting to be included in this thesis, since one of the objectives of 
this thesis is to provide a detailed survey of the properties of multi- 


variable root loci. 

We will discuss and prove three main results. First, it will be shown 
that the SISO root locus rule on the conservation of the sum of the closed- 
loop poles as the gain k is varied (see Section 2.5, rule (9)) generalizes 
to the multivariable case. This was first pointed out to us by Levy [27], 
and proved independently by Byrnes and Stevens [20]. Next, it will be shown 
how to compute the values of the gain k for which the root loci intersect 
the imaginary axis. This result is due to Shaked [28]. Finally, it will be 
shown how graphical bounds on the root loci may be constructed. This 
result is due to Owens and Field [29]. 

7.2 Results 

We now show that the SISO r. ot locus rule (Section 2.5, rule (9)) on 
the conservation of the sum of the closed-loop poles as the gain k is varied 
generalizes to the multivariable case. 


TTieorem 7.1 The sum of the closed-loop poles does not vary with the 
feedback gain k if and only if there are m infinite zeros, each of which 
has order not less than two. 


Relink . tis condition will hold if and only if all of the structure 
indices of G(s) at infinity are greater than or equal to two. 


We prove Theorem 7.1 by writing G(s) as an irreducible right matrix fraction 
description 

G(s) » N(s)D(s)“^. (7.1) 

The open*loop poles are then the solutions to 

DET (D(s)J » 0 (7.2) 

(see [23]), and the characteristic equation becomes 
A(g,s) * DET [gl - G(s)] = DET [I + kG(s)] 


DET [D(s) ^ kN(s)l . 
DET D(s) ^ 


(7.3) 


which may be written as 

m 

DET (D(s)l + kl^DET [0^ ... D^_^, Dj^^j ... Dj + ...» 0 

(7.4) 

where D^ is the i^^ column of D(s), and similarly for 

Now, DET (D(s)] is a monic polynomial of degree n, and the sum of its 
roots is minus the coefficient of s”“^. If and only if the rest of the 
terms of (7.4) are polynomials of degree not greater than n • 2, then the 
sum of the closed- loop poles does not vary with k. We now investi/iate the 
circumstances under which this is true. 

Let the column degre^ss of D(s) be [d^ ... d^], where d^ is the largest 
of the degrees of the elements of column If we assume that D(s) is 
column-reduced (see [12]), then we have 
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n ■ DEG [DET D(s)] « ♦ dj ♦ (7.S) 

Now, suppose that the column degrees of N(s) are [n^, ... n^]. If we have 

ttf i di - 2, i » 1, ... m (7.6) 

then the coefficients of all of the powers of k in (7.4) will be of degree 
not greater than n - 2. The condition (7.6) will be satisfied if and only 
if 


LIM sG(s) = 0. (7.7) 

S-H» 

Note, that the assumption of 0(s) being column-reduced is necessary for 
(7.6) and (7.7) to be equivalent. However, D(s) may always be made column- 
reduced by elementary column operations, and these may be included in the 
MFD (7.1). Both of these results are shown in [12]. 

Now let the Laurent series expansion of G(s) at infinity be 

G(s) » G. ♦ Gj + , . . (7.8) 

s * s^ ^ 

From Lemma 3.1 (see also [12], [23]), the number of first-order zeros at 
infinity is RANK Gj. If this is zero, is zero, (7.7) holds, and the result 
follows. 

Byrnes and Stevens [20] prove this result by using a different technique. 
The importance of this result is that it is a straightforward generalization 
of the corresponding SISO root locus rule, and as such it is worth noting. 

Let us now investigate the values of k for which the root loci inter- 
sect the imaginary axis. The importance of this for stability and gain 
margin should be quite apparent. 



Theorem 7.2 Let a closed-loop system be described by the state-space 
formulation (2.2), and let 


p rv^ ^ ^/princip&l minors of\ 
- it-^order i of kBKC - A^* 


i * 1 n. 


(7.9) 


Then if n is even the gains k at which root loci intersect the imaginary 
axis are solutions of the equation 


1 E2 E, 
0 1 E, 


E - E 
n-2 n 



- - 1 
2 ^ 


II ^3 ^5 


=n-l 0 


A similar determinantal equation applies if n is odd. 


Remark The polynomial equation (7.10) has degree 


n(n - 


(7.10) 


Theorem 7.2 plays the same role in the multivariable case that the 
Routh-Hurwitz criterion does in the SISO case, i.e. it can be used for finding 
the gains at which the system becomes unstable. In fact, (7.10) can be 
viewed as a multivariable generalization of a Routh-Hurwitz array. 

The proof of Theorem 7.2 is simple. We may write the characteristic 
equation in £ as 

A(k,s) « DET (si - A ♦ kBKC) « s" ♦ SE 4 Ch)s""^ - 0. (7.11) 

i«l ^ 

Letting s « ju and setting both the real and imaginary parts of (7.11) equal 


to zero gives 
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n 

I 


w ♦ ^(-1) E,. (k)w 
i-1 

n 



(k)a. 


n-2i«l 


0 


0 


(7.12a) 


(7.12b) 


if n is even and 

(n-1) 

2 

u” ♦ E (-1)4,. (k)«”’^^ - 0 (7.13a) 

i«l 

Injli 

E (-1)42. ,(k)u)""2^*^ » 0 (7.13b) 

i*l 

2 

if n is odd. Letting z « -u and using the resultant (see Appendix) to find 
values of k that will allow (7.12a.b) to have a simultaneous solution yields 
(7.10). 

Of course, we only need to search for positive real roots of (7.10), 
which simplifies matters considerably. In [28] Shaked noted that all positive 
real roots of (7.10) correspond to one of the following: the desired critical 
gains; the zeros of (7.10) that are symmetric with respect to the imaginary 
axis (a rare occurrence) ; and the loci passing through the origin. 

We now show how graphical bounds for the root loci may be computed. 

These may be useful if we desire only to have an approximate idea of where 
the loci are located (e.g. left half-plane versus right half-plane, for 
stability) . 


Theorem 7.3 Again using the state-space formulation (2.2), suppose 
there exists some real constants a, b such that 


4 

5 


aBKC bBKC 
•bBKC aBKC 


E(a,b) 


i. 0. 


(7.14) 


US 


Then the root locus lies entirely in the region 

R - {s: RE[s(a - jb)J < X(a,b)} (7. IS) 

where x(a,b) is the largest eigenvalue of the Hermitian matrix 

H(a,b) « i(a - jb)A ♦ l(a ♦ jb)A^. (7.16) 

Remark Note that (a,b) « (1>0) implies that the closed>loop system is 
stable for all gains if all of the eigenvalues of 1(A ♦ A*^) are 
negative. 

The following proof of Theorem 7.3 is borrowed from [29]. The charac- 
teristic equation (7.11) may be restated as 

(si - A ♦ kBKC)x(k) « 0 (7.17) 

U 

where x(k) is a non- zero n-vector with Hermitian transpose x (k), and where 
x^(k)x(k) is normalized to unity. Premultiplying (7.17) by x*^(k) yields 

s(k) - x“(k)Ax(k) - -kx”(k)BKCx(k) (7.18) 

and multiplying (7.18) by (a - jb), taking real parts, and writing 

x(k) » u ♦ jv (7.19) 


yields 


RE [(a - jb)s(k)l - xV)H(a,b)x(k) » -ktu”^ v''’]E(a,b)[“j. (7.20) 

Since H(a,b) is Hermitian and since E(a,b) ^ 0, we have 

RE ((a - jb)s(k)l < x“(k)H(a.b)x(k) < [x”H(a,b)x] 

■ X(a,b) (7.21) 


which proves Theorem 7.3. In [29] Owens and Field refine this result further, 
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giving simpler conditions for E(a,b) ^0, and discuss several applications. 
Plotting the bounds for several values of a and b can give a considerable 
amount of information on the whereabouts of the root loci. 

In [30] Brockett and Byrnes discuss the asymptotic behavior of multi- 
variable root loci from a geometric perspective. They consider the case where 
G(s) is non-square, and the case where the polynomial gain 

K(k) « ♦ Kjk ♦ ... ♦ K^k** (7.22) 

is applied to the closed-loop system in place of the gain kl. It so happens 

that if G(s) is non-degenerate , the asymptotic behavior is determined solely 

d 

by the highest-order term K^k , even if does not have full rank. G(s) is 
non-degenerate if and only if at least one closed-loop pole becomes infinite 
for any K(k) with d non- zero as k-^. If G(s) is pxm and has MacMillan degree 
n, non-degeneracy is generic if mp <. n. However, diagonal and block-diagonal 
transfer function matrices are degenerate. 

Brockett and Byrnes also show that if has full rank then the closed- 
loop poles that remain finite approach the open- loop finite zeros as k-^. 
However, both this result and the previous result depend heavily on the non- 
degeneracy of G(s) . For proofs and more details the reader is referred to 
[30]. 
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CHAPTER VIII 
Conclusion 


In this thesis the properties of multivariable root loci have been 
analyzed from a frequency-domain point of view. The behavior of the angles 
of arrival and departure has been sttJdied in considerable detail, and 
several methods of computing them have been presented. The problems of 
locating the break-in and break-out points and of characterizing the nuiaber 
of loci on the real axis, previously unexamined in the multivariable case, 
have also been addressed. Some methods for conq>uting angles and pivots for 
first-order asynqptotes have been given, and results for higher-order asymp- 
totes have been reviewed. At all times our objective has been to generalize 
the SISO root locus rules to the multivariable case, and as often as not 
such generalizations were found to exist. 

In Chapter III the angles of arrival and departure were analyzed using 
several different approaches. First, the state-space results of Shaked and 
Thompson were reviewed. Next, the general case of multiple poles and zeros 
was analyzed using the Smith-MacMillan form of G(s) , and equations for the 
angles derived. It was found that the angles are grouped into Butterworth 
patterns with orders given by the MacMillan indices of G(s) at the pole or 
zero in question, subject to certain conditions. Thus, in the most general 
case, the angles of arrival and departure are far more complicated for multi- 
variable systems than they are for SISO systems. 

We then specialized to the case of "simple** multiple poles and zeros, 
and found in this case that the SISO rules ^ generalize to the multivariable 
case. The reason for this is that non- simple poles and zeros have no 
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counterparts in the SISO case, so we can hardly expect SISO results to 
generalize to them. In fact, throughout this thesis, most of the failures 
of SISO root locus rules to generalize to the multivariable case can be 
traced to the presence of features not found in SISO systums (e.g. non- 
simple poles and zeros, branch points, multiple loci on the real axis, etc.)* 
Finally, the results of Sastry and Desoer on asymptotes, which involve solving 
generalized eigenvalue problems in Toeplitz matrices, were adapted to find 
angles of arrival and departure. 

The results of Chapter III were all conditioned on several generic 
assunqptions. However, with the exception of the first-order case, we were 
not able to relate these assumptions to each other and show that they are 
equivalent. More work needs to be done in this area, especially for the 
results derived from the Smith-MacMillan form. The cases where the assump- 
tions do hold should also be investigated. Reviewing the non-generic 
behavior of asymptotes (e.g. fractional orders; see [IS]) should be a 
useful guide in this endeavor. 

In Chapter IV the effect of branch points on loci on the real axis was 
examined. The "turnaround" behavior of loci at a branch point was considered 
in some detail. It was also shown that the SISO rules for break-in and 
break-out points generalize directly to the multivariable case. Some methods 
for computing both types of points, involving the resultant, were also 
presented. 

One feature of the siultivariablo root locus that definitely needs 
further investigation is the peculiar behavior that seems to be associated 
with branch points and branch cuts (see [S], p. 64). In this thesis we 
investigated only one small aspect of this issue. A characterization and a 


t 


I 


mm 




physical interpretation of this behavior would be of great help for under- 


standing and plotting the multivariable root locus, but achieving this 


will apparently be difficult. An interpretation of the branch points in 


terns of the state-space matrices (A, B, C) would also be useful. In 


contrast, there seeu to be little additional work that needs to be done 


on break points, except for finding easier ways of computing then. 


Multivariable root loci on the real axis were discussed in Chapter 


V. It was shown that, unlike in the SISO case, knowledge of poles and 


zero locations alone is insufficient for determining how many loci lie on 


the real axis at a given point. The reason for this is that in the multi- 


variable case several loci can lie on the real axis at a given point. 


whereas in the SISO case only one branch can lie on the real axis at a 


given point. An equation for the number of loci on the real axis at any 


point was derived. This equation requires only the computation of a few 


quantities involving the poles, zeros, and branch points on the real axis. 


A simplification of this equation would be helpful, but there are no evident 


approaches to take. There remains also the problem of extending this 


equation to the case of non-simple higher-order poles and zeros on the real 


axis, without bringing in the Smith-MacMillan results which, in our opinion. 


would make the general equation too cumbersome to be useful. 


The asyng)totic behavior of root loci was the subject of Chapter VI. 


Two methods of computing the angles of first-order asymptotes were given, 


as well as a method for computing the pivots of such asymptotes . It should 


again be noted that pivots in the multivariable case may be complex. 


Previous results on higher-order asymptotes were briefly reviewed, mostly 
for the sake of cosqileteness. In the context of the methodology of this 


I 

i 

) 

I 

thesis, these seems to be little edditionel work to be done here. Itowever, | 

i 

there ere many unanswered questions in the std)Ject of hi^er-order asymptotes, | 
particularly in cases where generic assumptions do not hold. Even the orders | 
of the asyi^otes are uncertain in som of these cases. 

In the SISO case the root locus is a handy tool for designing lead- lag 
compensators. It is well-known, for exaaq>le, that introducing a zero to the 
left of a branch tends to '*pull*' the branch toward the zero. A study of how 
these SISO tendencies generalize to the multivariable case would be very 
helpful. Of course, this is a very complex std>ject, but the results may 
prove to be worthwhile. Fo.v example, it might be possible to introduce 
a branch point and "turn around" an unstable asymptote. It is hoped that 
the results of this thesis will be helpful in studying the effects on the 
root locus of the Introduction of poles, zeros, and branch points in various 
locations. 
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APPENDIX 
The Resultant 


The resultant is a method for determining whether two polynomial 
equations have a common zero. It is constructed as follows. Let 


a(x) * ♦ . 

.. ♦ a_x" 
n 







(A.l) 

b(x] 

= b^ ♦ b.x ♦ . 
0 1 

.. ♦ b^x 

m 







(A.2) 

be two polynomials. Then a(x) and b(x) 

have a common non-constant factor 

if and only if 
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(A. 3) 

This result is proved in 

Walker [6] 

as 

follows. 

Suppose that a(x) and 

b(x) have a common non-constant factor (x - 

p)* 

Then we may write 


a(x) 

« c(x)(x - p) 








(A.4) 

b(x) 

« d(x) (X - p) 








(A.S) 


where 


• « • 


(A.6) 


c(x) * ♦ CjX + 




n-1 


<i{x) * + d.x ♦ ... + d ,x 

o i ni -1 


m-1 


(A. 7) 


have degrees one less than those of a(x) and b(x)> respectively. We may 
confine (A. 4) and (A.S) into 


d(x)a(x) - c(x)b(x) ■ 0 


(A.8) 


and substituting (A.l), (A. 2), (A.6), and (A. 7) into (A.8) and setting the 
coefficients equal to zero yields 


da - c b^ » 0 
o o 0 o 


da. + d,a„ - c b, - c,b„ * 0 
Ol 10 ol lo 


(A. 9a) 
(A. 9b) 


d_ ,a - c_ ,b = 0 
m-i n n-i III 


(A. 9c) 


which may be written as 


[do, di, ... '^2* *•• 


S ^1 ^2 
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0 . 0 a^ a, 
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"n 0 
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*>1 ^2 


b b. b^ ... 

m 0 1 z 
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n-1 n 
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i 0 
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^m ° 


(A. 10) 
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Equation (A. 10) will have a non*zero solution if and only if (A. 3) holds. 
This completes the proof. 

An inqportant application of the resultant is determining whether a 
polynomial equation has a multiple root. This issue arises in Chapter IV, 
since branch points and break points are multiple roots of 4(g,s) » 0. To 
determine whether a(x) « 0 has a multiple root, sicqily choose b(x) » ^ a(x). 
The left side of (A. 3) then becomes the discriminant of a(x). For quadratic 
a(x), 

2 

a(x) = ax + bx + c = 0 (A. 11) 


(A. 3) becomes 


DET 


c b a 
b 2a 0 
0 b 2a 


= -a(b^ - 4ac) = 0 


(A. 12) 


the familiar discriminant for the general quadratic equation. For the 
reduced cubic equation. 


a(x) * x^ ♦ px + q s 0 


(A. 13) 


(A. 3) becomes 
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DET 
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0 10 
p 0 1 

3 0 0 
0 3 0 
P 0 


27q^ + 4p^ « 0 


(A. 14) 
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which is also well-known. 

To aid in the computation of branch points for three-input-three 
output (m » 3) systems, the resultant for the general cubic equation 

3 2 

X +bx +cx+d*0 
is 

A s ISbcd - 4b^d + b^c^ - 4c^ - 27d^. 


(A. IS) 


(A. 16) 


